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ABSTRACT 


A software package was developed to perform conceptual 
design of submarines, using the Computervision CGP-=@2x 
Designer System, a turn-key computer aided design 
hardware and graphics software system. The philosophy 
behind the software package is to keep all major design 
decisions under the control of the design engineer, 
rather than embedding key decisions in the program 
algorithms. 


Modules are provided for calculating weight estimates, 


Principal characteristics and envelope geometry, 
resistance, weight and moment balance, and the 
equilibrium polygon. The package interfaces with a4 
pressure hull design module develaped separately im an 
O,.E. thesis by Marvin Meade. Interactive graphics are 


used where appropriate. 

The software requires a knowledgeable naval architect as 
the user, but does not require extensive knowledge of 
computers or computer aided design systems. 


Thesis Supervisor: Dr. David V. Burke 


Title: Professor of Ocean Engineering 
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INTRODUCTION 


GenmenAl BISCUSSION 


Computer Aided Design (CAD) has become an important 
tool in many industries over the past ten years. The 
capabilities of CAD systems are growing at a rate which 
defies efforts to catalog them. In the aerospace and 
automotive industries, as well as in large architectural 
firms, CAD has become the standard mode of design. 

In the area of naval ship design, CAD has been 
adopted more slowly, for a number of reasons. The 
comparatively low production rates bring the issue of 
cost effectiveness into greater question for ship design 
agencies. The tremendous complexity and scale of naval 
combatants makes the design of a flexible and 
comprehensive package of software a formidable task. 

Despite these problems, a number of packages for 
conceptual and preliminary design of surface combatants 
have been implemented in both government and private 


sector design organizations. The Packages used by 





gavernment agencies have tended to be  batch-oriented 
ship synthesis programs with Ligcles GY sno. araph ve 
output, and virtually no real-time interactive graphic 
design capability (the notable exception to this is the 
use Of interactive interior ship layout software). 

This lack of interactive graphics has been 
Partially hardware driven; the concern ror 
standardization and the lengthy process of government 
certification and procurement cause a substantial delay 
between the initial avallability and subsequent 
installation of new hardware. In addition, the CAD 
process is foreign to many of the prominent designers in 
the very tradition-oriented world of naval architecture. 
Many design procedures in naval architecture are 
somewhat subjective and difficult to quantify. 

A further impediment to the use of CAD has been the 
tendency of many early software packages ta usurp 
designer perogative by having key decisions in the 
design process "hard-wired" into the program algorithms. 
Such software does not utilize the expertise of the 
designer, and frequently precludes trading off design 
Parameters to achieve optimal designs. 

In the area of submarine design, the above problems 
are further accentuated. The volume of new submarine 
designs is even lower than that of surface combatants. 
Accuracy becomes critical} where feet and tons may “he 


acceptable units for tolerance in large surface ship 





designs, inches and pounds are frequently the desired 
accuracy level in submarine design. This is due to the 
hydrostatic requirement to balance both surfaced and 
submerged, as well as to the premium attached to 
internal volume and deck surface area. 

Submarine designers are fewer in number than their 
surface counterparts, and they often have many design 
relationships tied to subiective criteria CA Sas 
"designer’s eye") developed over years of experience. im 
program which does not allow the submarine designer to 
retain this subjective design flexibility is not likely 


to be accepted by the design community. 


SUBMARINE CAD? A HISTORICAL PERSPECTIVE 


Although Computer Aided Design has not been a 
Primary tool in submarine design, numerous programs have 
been developed to accomplish analytical chores that are 
subsets of the overall design process. These include 
Propellor design; resistance estimates, hydrostatic 
Calculations, and structural design. These programs have 
been primarily analytical in nature, providing noe 
interactive graphic capability, and are usually written 
for main-frame computer systems. 

Several packages have been written to achieve a 
complete iteration of the design process. These include 


the CODESUP program CijJ. It provides numerical, rather 





than graphic, output. CODESUR was developed by the 
Center for Naval Analuses in order to aid in projecting 
future submarine design characteristics for both our 
navy and its potential adversaries. This program was not 
intended to be a detailed tool for the naval architect, 
ands, unfortunately, hasS many key design parameters 
imbedded in the source code, which restricts its use in 
trading off alternative designs. Other packages "graft" 
together sections of existing designs, matching the 
largest hull diameter, to provide a composite design. 
ASSET C2j, developed bu Boeing Computer Services, while 
Primarily a surface ship design package, may be expanded 
to include a submarine design module in the future. 
ASSET provides some low-level graphics output, but is 
not a full, real-time interactive graphics package. 

THUS » there 1S a clear opportunity far the 
development of a submarine design package incorporating 
the use of the full capabilities of currently available 
CAD systems. This thesis represents an exploration of 


such a design package. 


Pel LOS@En Yer ORMIIHE THESIS PROGRAM 


The software developed in this thesis accomplishes 
ae e Pass” through the submarine design process at 
the conceptual design level. In Congtierc ron with the 


Gonewrcrent Ocean Engineer Pies Ss written by Marvin 


Meade, the package includes calculation of weight 
estimates, envelope and pressure hull geometry, 
Pesi=tance (speed and power), basic hydrostatics, weight 
and moment Balance, and determination of the equilibrium 
Se Log @t) 

tae primary philosophy of this package is to 
develop it for a user who is a competent naval 
architect, but who may not necessarily be an expert in 
the use of computers. Crucial design decisions are left 
to the user, rather than being embedded in the program 
code. Opportunities are provided to override program 
algorithms where necessary to implement the desires of 
the individual designer. In addition, if the designer 
does posess some knowledge of CAD systems or computer 
programs, it 1S easy to exit the package at appropriate 
Points and tailor the design to provide more detail or 
alternate geometry and analytical procedures. 
Interactive graphics are used to provide a clear 
Visualization of the design under development. All 
program input parameters are couched in standard naval 
architecture terms to facilitate the user’s interaction 
With the package. The title chosen for the software 
Package is CADSUP (Computer Aided Design of Submarines). 

Input to each program module is explained in the 
chapter documenting the particular module. Samples of 
Program output are provided in Appendix I to this 


thesis. 
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HARDWARE : 


THE COMPUTERVISION DESIGNER SYSTEM 





GENERAL DISCUSSION 


The optimum choice of hardware for implementing the 
Philosophy of this thesis was one of the several 
high-quality "turnkey" CAD systems, all of which permit 
the use of interactive graphic design techniques. 
Graphics packages for main-frame systems are available 
from several prominent suppliers, but are a compromise 
at best, since the processors they run on are not 
optimized for graphics-intensive applications. 

Suitable systems are marketed by several suppliers, 
Pa vein Computervision, Applicon, Bendix, 
McDonnell—-Douglas Automation Division, and others. For 
this thesis, Computervision Corporation, of Bedford, 
Massachusetts, offered free system time, instruction, 


anGeoartnwenm= technical aid to develop the program. 


THE COMPUTERVISTION DESIGNER SYSTEM 


The system which Computervision provided time and 
instruction on is their basic "Designer System", elt 
around the CGP-=@OxX processor. The heart of the sustem 
consists of the processor, a high-speed tape drive, and 
one or two 300 megabyte hard disk units. The CPU is 
available with several increments of random access 
memory (RAM). The particular system used for this thesis 
was configured for approximately 1.3 megabytes of core 
memory. The architecture of the CGP-2@@xX is 16-bit. 

The Designer System utilizes two distinct operating 
sustems. Basic file management and syustem—-level 
housekeeping chores are handled at the "OS" level, which 
is Similar in use to the operating systems on typical 
MmMini-~computer installations. For the graphics operating 
system, Computervision uses the "CADDS 4-X" environment, 
which 1S tailored for graphics, and automatically 
interfaces with a Powerful data base management system. 
Graphics commands are entered in a simple verb~-noun 
Suecesesticn as “INSERT LINE” or “INSERT SPLINE”. 

Input/output (1/0) is handled by several devices. 
For system administration and text file input, small 
alpha-numeric terminals are provided. For graphics [/7/0, 
high-resolution color or monochrome terminals are 


Provided as part of a comprehensive workstation, which 





includes a thermal printer for WOrking quality text and 
graphics printout. Also provided at each workstation are 
a digitizing tablet and a display control device. This 
device controls background and graphics intensities, the 
number of lines of text displayed (4 or 24 lines), and 
waattercOntrol of the display$ zooming; rotation, and 
scrolling. Sutput of finished quality is available from 
line printers, color pen plotters, and black and white 


electrostatic plotters. Large format output is supported 


by the plotters. 


PROGRAMMING LANGUAGES 


———es eee 


Several languages are provided for use with the 
Designer System. The primary iaeauase 1s (FPortran—os a 
subset of ANSI Fortran. System subroutines are provided 
for interacting with the graphics and data base 
management operating systems. Programs entered in 
Fortran-S execute with the greatest speed, and allaw 
creation of new graphics commands. Because of limited 
CPU space, the utilities required for this level af 
programming, which are often provided on main-frame 
sustems, must be separately programmed for each module. 
Such utilities include linking and loading capabilities. 
Pecause of the length of training required to do 


competent programming at this level, this thesis was 


programmed in an alternate language. 





Two macro languages are provided for the system. 
Poetic bow Greation Of execute files using standard 
graphics commands. VARPROZ is a non-compilable macro 
language which Supports Fortran-like 70 and 
computational statements. It is a flexible and effective 
language, but executes rather slowly. NEWVAR, the other 
macro language provided; is compilable, and therefore 
executes far more rapidly. The I/79 for NEWVAR is 
cumbersome and difficult to format. All I/O must be 
handled in text string form, and converted within the 
program to numerical values. Despite this shortcoming, 
the execution speed of NEWVAR led ta its adoption (“far 
use in this thesis. Ultimately, “the algorithms in this 
thesis should be translated to Fortran-S for maximum 


speed and [70 flexibility. 





PRELIMINARY WEIGHT ESTIMATES 





IMPORTANCE OF WEIGHT ESTIMATES 


The algorithms employed in this thesis assume a 
weight-driven submarine design. This assumption means 
that the envelope size calculated for the submarine will 
contain enough valume to displace the same amount = of 
water, when submerged, as the weight of the envelope 
displacement for the submarine. If the desired design 1s 
known ar suspected to be volume-limiteds, an arbitrary 
excess of lead may be specified during the weight 
Calculations, thus driving the program to calculate a 
larger envelope size. The designer may translate this 
additional weight into volume manually to allow for 4a 
volume-limited design. 

Since the estimated weight data will drive the 
basic dimensions of the submarine, the weight estimating 
algorithms must be made as accurate as possible. The 
designer must carefully consider the input data to the 


weight module, and, where appropriates may decide 16 


override the programmed algorithms and specify the value 
for a particular weight group. The program 1s designed 
to automate the calculations, and not the design 
Judgement, for the user. A knowledgeable user may easily 
modify the algorithms themselves to suit particular 
needs. 


WEIGHT ESTIMATING ALGORITHMS 


The accounting system chosen for weight estimaticn 
is the SWEPS (Ship Work PBreakdown Structure) canvention 
used by the Naval Sea Systems Command for ship weight 
records C3]. This system consists of seven major 
numbered groups for all ships. For submarine weight 
reports, several special weight categories are added. A&A 


brief stummary of the weight groups is as follows: 


GROUP DESCRIPTION 
I Structural 
Te Propulsion 
ei Electpiceat 
IV Combat Systems 
V Auxiliary Systems 
VI Qutfit and Furnishings 
VII Armament 
Al Sint Gipmeeo vy al 
Lead Total lead 
- Ai + Lead 


Variable Load (fuel, stores, etc.) 


Sos 


NSC Normal Surface Condition 
Reserve Buoyancy Main Ballast 

Submerged Disp. self explanatory 

Free Flaad Non-buoyant flooded volume 
Envelope Disp. Submerged Disp. + Free flood 


The algorithms chosen for the weight estimating 
module take standard design parameters as input, where 
possible, and provide an output in terms af a percentage 
of Al weights, NSC, or submerged displacement, as 
appropriate. These percentages are then combined with 
the discrete weights that are calculated or input for 
selected weight groups to determine the Al weight total, 
lead, NSC, and other weight parameters. The designer has 
the option of specifying weights far groups II, IV, VI, 
and VII, which override the programmed algorithms. The 
following section discusses the chosen algorithms for 
Calculating each weight group. Throughout this thesis, 
the following conventions will be used to indicate 


mathmatical operations: 


— = acielmt ian 

"—" = subtraction 
eee eli ltipliacation 
"/" = division 

"A" = exponentiation 


GRewWP 1: 


Weight group I is determined as a percentage of 


NSC. In general, 
ZW1i=f( depth, material, size). 

The size dependence is accomodated by Calculating 
Wi as a percentage of NSC. The algorithm implemented in 
the program is: 

At Ole ee DEPTH. 

Cl and C# are dependent on hull material used, and 
are determined parametrically from data on Past 
submarine designs [4,5]. The DEPTH parameter is a user 
input of maximum Operating depth in feet. The accepted 
design practice of designing for a collapse depth of 
150% of operating depth is implemented through the 


choice of Cl and Cz. 


GROUP II: 

Weight group II is calculated directly, in tons; 
rather than as a percentage. The number of parameters 
involved in calculating WE is greater than for any other 
weight group. 

We=fChorsepower, propulsion type, battery type & 
capacity (capacity=endurance) J] 

The horsepower variable for a nuclear plant 15 
SimPly shaft horsepower, while, for a non-nuclear Plant, 
two separate horsepowers are needed; shaft horsepower 
(usually the horsepower of the electric final drive 
motor) and charging horsepower (power of the prime mover 


used to drive the alternators or generators oe provide 





electrical power to recharge the storage batteries). The 
weight of the charging prime mover is a strong function 
Of technology chosen; ie diesel, wankel, etc. 

The calculation of storage battery weight for a 
non-nuclear plant is dependent on battery tupe, Capacity 
mn Ellowatt—-hours at a high discharge rate ("sprint"), 
and capacity in WKHW-H at a low discharge rate (endurance 
Power). The options provided for battery type in this 
program are nickel-cadmium, improved lead-acid (German 
VARTA type), and "standard" Lead aeac ome CU .S. Trident 
type). Constants determined from the energy densities of 
these battery types at the two discharge rates are used 
to determine battery weight (the larger of sprint or 
endurance rate battery weights is chosen as the final 
battery weight). 

Electric motor weight 15 determined from a 
multiplier based on a modern, air-cooled DC mator. The 
OoOptian exists to provide other multipliers based on 
liquid cooling, superconducting technology, etc. 

For the non-nuclear plant; an additional option 
exists to provide a lumped adjustment to weight group I[] 
to provide for unconventional propulsion technology 
(fuel cells, Stirling, low-power nuclear, SiGe 7. This 
final adjustment allows the designer maximum flexibility 
ase geoup [iowithout making modifications ™to the 
program code. 


For a "straight" nuclear plant, 


We=(CL¥SHP) /ED log (SHP) 14C2 

where SHP=shaft horsepower 
Ci is from parametric data. 

For a non-nuclear plant, 

WH=WE+WM+WER+dWe 

where WE=charging prime mover weight 
WM=propulsion motor weight 
WR=battery weight 
dWe=adjustment to WE (user input) 

HE=C5* CHP 

where CHP=charging horsepower 
C3 is dependent on prime mover type. 

Wh=C4*SHP. 

WR=Maxl (CS*KWHS) » (C6E*KWHE ) J 

where KHWHS=capacity for "sprint" 
WWHE=capacity for endurance 


C5,C6 are dependent an battery tupe. 


GROUP TIT: 

Although weight group III would intuitively appear 
to be a direct function of installed generating 
Capacity, a study of historical design data [51] reveals 
that the group III weights are, in fact, very closely 
approximated by a straight percentage of Al weights. 
This anomaly is best explained by the tlarge amount of 
this weight group attributable to the electrical 


distribution system spread throughout the submarine. 


This power distribution system weight is a direct 
Function of the size of the submarine, since the amounts 
of cable, cannectors, and distribution panels vary with 
size of the vessel. Weight group III, consequently, is 


calculated as a percentage of Al weights: 


GROUP IV: 

Determination of group IV weights is very difficult 
to implement in an "automatic" mode. The variety of 
combat system components and great variation in 
equipment weights makes the selection of an algorithm 
extremely difficult, as most schemes provide accurate 
results for a limited number of combat suites, at best. 
It is therefore highly recommended that the designer 
specify a group IV weight based on his off-line 
determination of the desired combat systems and 
associated weights. To permit a rough approximation 15 
acceptable for an early design iteration, the program 
contains an algorithm Based on historical data as a 
percentage of Al weights for either fast attack or 
ballistic missile submarines: 

AZW4=Ci i=1,4 
where Ci is based on submarine type (55, SSP, 


SSN, or SSBN). 


GROUP V: 

Weight group V is accurately determined as a 
function of submarine size, since the auxiliary systems 
required for ae vessel are directly proportional to 
displacement and internal volume. Group V is thus 
calculated as a percentage of Al weights: 

%W5=Ci t= lene 
where Ci 1S based on submarine type (attack or 


ballistic missile). 


GROUP VI: 

Group VI weights are a function of crew size, and, 
to some extent, typical mission duration. The multiplier 
for ogrompe VI weight determination is based on 
habitability standards commensurate with the last two 
classes of attack and ballistic missile submarines, 


respectively. The weights are directly Seer] wee an 


tons: 

W6=C1*NP 

where NP=number of personnel in the crew. 
GROUP VIT: 


The estimation of armament weights 1s nearly as 
Gittict@ee as the task™of determining group IV weights, 
Since the possiteie combinations of existing and future 
weapons and associated launchers are infinite. Again, it 


is highly recommended that the designer specify a value 





oe 


° . ——_ 
Bee ee 


for weight group VII based on his off-line determination 
of the weapons payload. As was provided in the case of 
combat systems, a rough estimate of group VII 1s 
available from a percentage of Al weights. The 
percentage constant was determined from historical data, 
and separate values are provided for attack and 


Ballistic missile submarines. 


AW/=C1. 


OTHER WEIGHTS: 

The amounts of lead, variable load, reserve 
buoyancy, and free flocd are provided as fractions by 
the designer. The following ranges are appropriate: 

Lead CoO 7—all) «Al 

Variable Load (.04-.@7)*NSC (nuclear) 
(OSS. Teo A nOna= Mitel ear) 

Reserve Buoyancy Gries eis eS & 

Free Flood (2a el) e CSolbmner ged Dicei. 

The operator, however, may specify his own values 


as desired. 


CALCULATION OF WEIGHT SUMMARY 


When either percentages or actual values have been 
Calculated or specified for each of the weight Groups, 
the calculations for the weight summary are processed. 


At a minimum, actual weight values etal tons have been 





determined for weight groups II and VI. Additionally; 
the designer May have specified weights for weight 
groups IV and VII. As an example of the algorithms used 
for final weight calculations, assume that only W2 and 
W6 have been determined in tons. The other permutations 
of final weight calculations are Similar, but allow for 
dicrete weights, vice percentages, for either W4, W7, or 
both. For the case where only We and W6& are knowns the 
sequence is as follows: 

(All "ZZ" figures used are converted to fractions; 
AX=7XK%/100) 

Ki=(ZW1i+ZW1#27ZLD)/01-ZVL)+2WS+2W5 

where “%ZLD=lead fraction 

%VL=variable load fraction. 

AL=(WE4+W6)/C1-K1-2ZW4—-ZAW7) 

WS=LZWS*A1 

W4=7W4#A 1 

W5=Z2w5*A1 

W7=2W7#A1 

Lead=/LD*A1 

A=A1i+Lead 

NSC=A/(1-ZVL) 

Variable Load=“VL*NSC 

Main Ballast=ZRPY#NSC 

where ZRBY=fraction of reserve buoyancy 

Submerged Disp.=NSC+Main Ballast 


Free Flood=ZFF*Submerged Disp. 





where Z“ZAFF=free flood fraction 
Envelope Disp.=Submerged Disp.+Free Flood. 

The calculated weight summary is displayed on the 
screen with a summary of user inputs, at which time the 
designer may accept or reject the results. ii ene 
results are rejected, the program loops back to the 
input section for the weight module and another 
iteration is started. If the results are accepted, they 
are written into an output file for record purposes; and 
various input and calculated parameters are written into 
several intermediate files for use in later program 
modules as the design is completed. 

Because of the high degree of importance attached 
to the calculated weight data, future refinement of 
Program algorithms should include a priority effort to 
accomplish two major tasks: 

(1) Better accuracy of calculated weights. 
(2) Further subdivision of weights (at least 


partially to the three-digit SWES level). 
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PRINCIPAL CHARACTERISTICS 


AND ENVELOPE GEOMETRY 


PoC LRA CHARACTERISTICS 


The results from the weight module pravide the 
first important input to the determination of the 
remainder of principal characteristics, which are 
geometric in nature. Since this thesis implements a 
weight-driven submarine design process, the next design 
task is to create an envelope of the appropriate volume 
to contain the envelope displacement: 

Envelope Volume=35*Envelope Displacement. 

The product envelope defined by this module of the 
program will exhibit the following geometric principal 
characteristics (abbreviations and symbols in 
Parentheses are program variable names): 

Length Overall (LOA) 
Length of Entrance (Lf) 
Length of Run (La) 

Mid Pody Length (Lm) 


Diameter (D) 





Length/Diameter Ratio (I./D) 
Prismatic Coefficient (Cp) 
Bare Hull Surface Area (SF). 

In addition, once envelope geometry is defined, the 
Following calculated hydrostatic characteristics may he 
found: 

scubmerged Longitudinal Center of Puouancy 
(EGERSUR.) 

Draft @ Normal Surface Displacement (TNSC) 

Longitudinal Center of Buiyancy @ NSC (LCBNSC). 

These parameters, then, when coupled with the input 
parameters delineated for the weight estimation module, 
constitute a set of Principal characteristics for the 


design. 


Dee RMONG Shae ENVELOPE 


A number of algorithms have been suggested to 
calculate the geometry of the envelope for submarines. 
These range from a "Cut-and-paste" approach using past 
designs C2 J co very sophisticated Polynomial 
determination of offsets [6]. In between these extremes 
are a number of methods utilizing simple thumbrules and 
Parametric relationships among CPs, D, a7 Wr and 
displacement €4,7 1]. 

Both the cut-and-paste and thumbrule approaches are 


somewhat imprecise in allowing the designer to Fate cs 


Gontrol the geometry of the resulting hull Pee Ms The 
polynomial methods, while very precisely controlled, 


require the designer to specify design parameters which 


Seeemea. COUGhKGG in, or clearly related ies traditional 
haval architecture vocabulary. The most desirable 
algorithkim, then, nie allow strict eoncrod ot the 


resulting geometry while accepting input Parameters 


which are common naval architectural terms. 


CHOSEN ENVELOPE GEOMETRY ALGORITHM 


In choosing the algorithm implemented in this 
thesis, two restrictions were accepted for the sake of 
Simplicity. The envelope shape is confined ta circular 
cross-sections, and an analytically smooth envelope 
cutline is calculated. Figure 4.1 demonstrates the 
nomenclature and subdivision of the envelope calculated 
with the chosen algorithm. 

The designer 15 prompted ior the following 
geometric input parameters: 

Length of Entrance 
Length of Rttn 
Forebody Prismatic Coefficient (CpF) 
Aft Body Prismatic Coefficient (CpA) 
Maximum Diameter (D). 
The algorithm then defines an envelope consisting 


of a hemi-~elliptical forebody, a eu a nGic leat midbodu, 





and a Parabolic afterbody. The chosen algorithm 1S a 


Moegification of one suggested by Jackson [£41]. 





a i nt ee nr ee ne C) 








Figure 4.1: Nomenclature for determination of envelope 


geometry. 


os 
= LD = 


The coefficients Nf and Na are essentially fullness 
coefficients for the fore and aft bodies, respectively. 
ineewerder §toO facilitate input in traditional naval 
architectural terms, these fullness coefficients are 
calculated from the values specified by the designer for 
prismatic coefficients of the fore and aft body sections 
of the envelope: 

NF= 71.0477*CpF45 +65.8107*CpF44 —-483.3037*CprF43 
Zooey er = Cpr oS —281.7224*eCprF 449.5876 
Where Nf is the forward fullness coefficient. 
Na= 379.6546*CpA%5 -938.4708*CpAr4 
+944, 9853*%CpASS -471.Q0872*CpA%S +119.1465*CpA -11.3454 
Where Na is the aft fullness coefficient. 

These polynomial relationships are valid for values 
Pactra. eo-., O72, and for Cpa from .3333-.7111. 
These ranges substantially exceed those for actual 
submarine designs. The relationships were derived by 
using known (hand-verified) equivalencies of Prismatic 
and fullness coefficients to set up matrix equations 
with the polynomial coefficients as the unknown column 
vector. Gaussian elimination was then used to solve far 
the unkown coefficients. Additional hand-calculated 
equivalencies were then utilized to verify that, within 
the specified ranges, the polynomial relationships yield 
results correct to the fourth decimal place. 

The "X" coordinate for the overall envelope runs 


from @ at the forward perpendicular to LOA at the aft 





perpendicular. In the formula hor the forebody 
hemi-ellipse, a local X coordinate (Xf) runs from @ at 
the end of the entrance length (Lf) to Lf at the forward 
perpendicular: 

Xf= Lf-X 

In the aft body parabolic section, a local xX 
coordinate is again used: 

Xa= X-(LF+Lm) 

Note that Lm, the mid body length, is a calculated 
Parameter, and not an input. The calculation of Lm will 
be described later in this chapter. 

The hull radius, r(X), for the elliptical forebody 
is calculated as follows: 

eer REET -CXE/L TANF ACI/SNE > J 
where R= D/z 

The hull radius within the mid body is simply: 
r(X)= R 

The full radius within the aft body parabolic 
section 1s: 

r(X)= R*C1-(Xa/La) “Nal 

The offset r(X) 18s calculated for each of £1 
stations evenly spaced along the hull length. Obviously, 
r(X) at both the forward and aft perpendiculars 1S zero. 
When the offsets for the #1 stations have been 
determined, the sectional area of each station is 
Calculated: 


SACX)= Pi¥r(X) 42 
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The station sectional areas are then integrated 
along the length to determine the volume enclosed by the 


envelope: 
LOA 
VOL= / SACX) dX 
v 


frtsweand abieother integrations in this thesis are 
done numerically, using Simpson’s Rule. 

Qn the first pass through this process, the mid 
body length 1S zero, and LOA= Lft+tLa. The resulting 
volume, then, is the total volume for the fore and aft 
body sections alone. To determine the proper mid body 
length, the volume determined above is subtracted from 
the required envelope volume, and Lm is calculated as 
the length of a cylinder of radius R containing the 
remaining volume: 

DV= Envelope Volume-VOL 
m= OP 1 eR A) 

If DV is negative, an error message is printed, and 
the designer is prompted for new geometric input values. 

After Lm is calculated, the final value of LOA is 
determined: 

LOA= Lf+Lmt+La 

ff ecoeeionamearcas for the resulting hull form are 
Calculated at the new £1 stations; after which 
circumferences and bare Full surface area are 
determined: 


Gueecx)= 2¥Pisrt x) 


“; 
= .* ee 
a 


LOA 
SF = 7 Caretx) dx 
Q 


The overall prismatic coefficient is calculated by 
comparing the envelope volume with the volume of a 
cylinder (Vcyl) of length LOA: 

wee Pie R4aelLlOA 
CPp= Envelope Volume/Vcyl 
These calculated parameters, when taken together 


wi th the Tmt values, constitute the geometric 


definition of the envelope. 
CALCULATION OF PASIC HYDROSTATICS 


In order to later balance the submarine, it els 
important to accurately determine the longitudinal 
centers of buoyancy in the submerged (LCPBSUB) and normal 
SCliriaee wGLCRNSC) comditiens. Calculatvem of LCRNSC 
requires the determination of the draft at NSC. LCR’s in 
both conditions are calculated by integrating the first 
longitudinal moment of the sectional areas. 

Using the sectional areas for the fully-submerged 
hull as calculated above, LCBSUP is determined: 

LOA 
LCRSUPR= [ X¥SA(X) dx /VOLENV 
@ 
There are several Possible approaches to 


determining the draft at NSC. When the calculations are 


done manually: it is customary to choose at least three 
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points, commonly at drafts equal to D, D/#, and D+R/2, 
and determine the volumes and resulting displacements 
for each of these drafts. A curve of displacement versus 
draft 1s then constructed, after which the cleat 
corresponding to the NSC displacement is read from the 
curve. While this method is expeditious and acceptable 
for hand calculations, relatively large errors may 
result when the designer is not experienced enough to 
draw the proper curve through the three points. Use of 
additional points improves the accuracy, but the manual 
calculations are tedious and time-consuming. 

With the speed and accuracy of the computer 
available, a much more precise determination of draft at 
NSC ¢(TNSC) is possible. A potentially important factor 
that is often overlooked in manual methods is the volume 
of the free flood, which 1S contained in the envelope, 
Tanase ae Sotoniftrigant portion (:'-104) of the envelope 
volume. Study of existing designs indicates that 
approximately 8@% of the total free flood volume is 
submerged at HNSCe Consequently, the "target" 
displacement in determining TNSC is not simply NSC, ee lite 
rather includes this portion of the free flood: 

NSCVOL= 35*(NSC+.8*Free Flood) 
(for=wour free flood). Where 
NSCVOL is the desired volume at TNSC. 
The method chosen to determine NSC draft 1s a 


bracket and halve iteration on the variable TNSC, 





beginning at TNSC= 3*R/2. The allowable error in NSC 
volume 1s specified as @.5% in the current program. This 
tolerance may be decreased, with additional iterations 
required for convergence. The value of @.5% was chosen 
as sufficient for the conceptual design phase while 
minimizing the number of iterations required. The 
algorithm employed is as follows: 

DV= .@@0@5*NSCVOL 

NPLUS= NSCVOL+DV 

NMINUS= NSCVOL-DYV 

DR= .5*R 

mNSC=R-EDR 

For each trial TNSC, the value of each station 
radius, r(X), is compared to DR: 

Tf (DR)er(X), then r(X) and SA(X) remain 
unchanged. 
ijeecpkR yar (Ad, then. 
SA(X)= SACX) -Er(X)42 *ARCOS(DR/r(X)) —-DR 
7 oa eee DR“ S) SL 34 
(Zero trim is assumed). 
LOA 
VOL= [save cIX 
@ 

If YOLSNMINUS and VOL<NPLUS, then the correct TNSC 
is determined. If VOL“NMINUS, then DR is increased Oy 
_S*DR. Similarly, if VOL>NPLUS, then DR is decreased by 
the same increment. In either case, thes ad gdoarithm 25 


repeated until the value of VOL is acceptable. When the 


correct TNSC is found, LCBNSC ic calculated: 
LOA 
ECan. — / X*¥SACXK) dX /SNSCVOL 
Q 

The calculated results are displayed on the screen 
For the designer’s approval. If the results are deemed 
unsatisfactory, the designer may reject the envelope 
geometry and return to the input section for another 
envelope design iteration. 

If the results are accepted, all of the specified 
input values and calculated data are written to an 
output file for record purposes, and appropriate data 


for other program modules are written to Passing files 


for later use. 
mem G lime ENVELOPE 


An accurate Plot of the envelope iS imperative, 
Since the outline displayed on the monitor will guide 
the designer in his interactive design of the pressure 
nok 1. The 21 stations determined thus far are 
insufficient to fit an accurate curve, since the 
arbitrary equally-spaced stations may not fall at the 
Points of greatest hull curvature, or at the Junctions 
of the mid body with the fore and aft body sections. 

To adequately portray hull curvature and 
transitions at section junctions, the program calculates 


13=Z uneauallu-spaced offsets to be used in plotting the 
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envelope. Offsets are most closely spaced near the 
forward perpendicular and at the transition points of 
the hull. These points are then used to fit a Bezier 
sPline curve which forms the profile outline of the 
envelope. The plot is automatically scaled to make full 
use of available screen = area; and both vertical and 
horizontal scales are drawn as a reference for the 
designer. 

The two-dimensional envelope outline is rotated to 
form a surface of revolution, and a mesh pattern is 
superimposed on the surface to aid visualization of the 
resulting three-dimensional hull form. In addition, an 
isometric view is displayed below the profile far 
perspective. 

The resulting display of two views of the outer 
hull constitutes the graphic output of the envelope 
geometry module. The designer is now ready to proceed to 


the speed and power calculations. 


SPEED AND POWER 


SENERAL DISCUSSION 


The calculation of resistance for the envelope 
geometry developed in the previous module 1s a 
straightforward analytic process. The formulae employed 
are the familiar drag estimation relationships from 
basic hydrodynamics. 4, proposal heir future 
implementation would be the addition of ae propellor 
design module to accurately predict the propulsive 
coefficient, in conyunction with the envelope geometry. 

Lo facilitate endurance calculations Con 
non-nuclear designs, the calculation of total electrical 
load at each speed is included. This total load figure 
must include all non-propulsion loads (hotel and combat 
system) as well as those required to propel the 
submarine. 

With the data passed from the envelope geometry 
module, an estimate of appendage drag may be made. This; 


when combined with hull form and surface area data, 
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completes the input to the speed and power alaori cht. 


ESTIMATION OF APPENDAGE DRAG 


Study of existing submarines reveals that appendage 
drag comprises thirty percent or more of the total drag 
on the submarine for most designs. In the later stages 
of the design process, detailed consideration should be 
given to the configuration and size of each appendage to 
improve the accuracy of resistance estimates. At the 
stage of conceptual design, however, details concerning 
appendages are often uncertain. Many designs are 
formulated without final selection of control surfaces; 
for example. The required shape, location and sizes of 
surfaces for cruciform, Xs, and inverted Y sterns may 
exhibit some variance, and the final configuration may 
not be chosen until much later in the design process. 

In view of this uncertainty, a reliable rough 
estimate of total appendage drag is desirable for the 
conceptual design. This thesis incorporates the method 
suggested by Bukalov [8]. Pukalov studied numerous 
actual submarine designs; and developed an algorithm 
based on this parametric study. This algorithm may he 
reduced to a single formula dependent on length overall 
(LOA) and maximum diameter (D). The appendage drag (Da) 
is then calculated as follows: 


Ba— |. O9065E—-3*L0A*Dt+11.225 


CHOSEN ALGORITHM 


The basic relationships implemented in this program 
module may be found in numerous introductory texts 
C9,1@]. The interactive input from the terminal consists 
of the following items: 

Propulsive Coefficient (PC) 
eee eet seon Loads in WW (CHIL) 

The remainder of the required input is passed from 
the envelope geometry module, and consists of: 

Prismatic Coefficient (Cp) 
Length Overall (LOA) 
Maximum Diameter (D) 

Bare Hull Surface Area (SF) 

The speed and power data 15 calculated and stared 
in an array indexed on speed in knots (VE) from zero to 
forty knots. After appendage drag is estimated in the 
manner previously described, the calculations proceed in 
an iterative loop: 

Re(Vk)= (VE#1.689*LO0A)/1.27908E-5 

Where Re is the Reynolds Number. 

Cia a/ lt Log (Re) —-2 aes 

Where Cf is the frictional drag coefficient, 
calculated by the ITTC convention. 

@e—vefel 1.5%(D/LOA)*%1.5+7*#(D/LOA)%3] 


tog ( Cp—. 6) 
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Where Cr is the residual, or form, drag 
coefficient, calculated after Jackson [£41], and the 
last term 15 a correction for submarines with 
considerable parallel mid body. 
Ge=—Cr+Cr+. 000.5 
Where Ct is the total drag coefficient for the 
bare hull, and the final term is the correction 
Factor from tow tank studies. 
HP(VED= CI/PC)*. GOB 72e¥VE 43 ( CtXSF+Da) 
Where HP(VkE) 15 the drag in horsepower for the 
given speed; Vk. 
KWOVED= . 745 7*HP CVE) +HL 
Where KWC(VE) is the total load in kilowatts, 
including non-propulsion loads; far the speed, Vk. 
These calculated results for drag in horsepower and 
total load in kilowatts are displayed on the terminal 
screen for designer inspection. From these results, the 
designer may determine the maximum speed of the 
submarine for the installed shaft horsepower, as well as 
the submerged endurance at any speed for installed 
energy storage capacity. If desired, the module may he 
re-run for various values of propulsive coefficient and 
non-propulsion loads, until a satisfactory result is 
obtained. 

When the designer accepts the calculated values, 
The results are written into an output file fie 


hard-copy record purposes. 
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PRELIMINARY BALANCE 





GENERAL DISCUSSION 


When the designer has completed the speed and power 
module, he will proceed to the design of the pressure 
hull, Which 1S aA process covered separately in a 
concurrent Ocean Engineer thesis by Marvin Meade (Cild. 
With a graphic display of the envelope and pressure hull 
to serve aS a visual aid, the designer will enter the 
subject module for this chapter toa accomplish a 
longitudinal weight and Buoyancy balance. 

The purpose of the preliminary balance is to 
determine the locations of lead and main ballast tanks, 
in order to confirm the feasibility of the pressure hull 
design. The required lead (Pb) longitudinal center of 
gravity (LCG) is achieved by balancing the longitudinal 
moments of all weights except the main ballast tanks to 
create am attitude of zero trim in the normal surface 
Eencdienoan <NSC). 


Following the determination of lead placement, an 


~ 


Se 


additional moment balance is computed to locate the main 
ballast tank LCG in order to achieve zero trim in the 
submerged condition. 

inpeseGiecton tO the longitudinal balance, a vertical 
weight balance 1S computed to provide an estimate of 
transverse Stability while Submerged. This estimate is 
in the form of BG, the distance from the vertical center 
Of buoyancy to the vertical center of gravity of the 
Submarine. This parameter should have a value of about 
one foot for an optimum design. To simplify 
calculations; no trim 1S allowed in the surfaced 


condition. 


CALCULATION PROCESS 
The first set of input data required for balancing 
the submarine is passed from previously completed 
modules. The data which 1s read in from earlier 
calculations is as follows: 
41116 weight categories from the Weight Module 
(Wi-W16) 
Longitudinal Centers of Buoyancy at NSC and 
submerged (LCBNSC and LCPSUB) 
Maximum Envelope Radius (R) 
Length Overall (LOA). 
The first set of keyboard input is the LCG and VCG 


information for the seven primary weight groups and the 


Variable load: 
mee 1) —1CG(7) 
WeGti=/CG( 7) 
ECGCVt) 
We ( Vs) 

Although detailed arrangements have not, as yet, 
been addressed in the design process, the designer must 
have given enough thought to rough arrangements ta 
Provide these location data Lor the aggregate weight 
subdivisions. If, when detailed arrangements Are 
completed at some later time, weight group LCG and VCG 
locations are signifigantly different from the input 
provided for this module, the balance process must be 
repeated. 

The program now calculates arrays of longitudinal 
and vertical moment data: 

LM(i)d= LCG(1i)#W(1) 
VM(1i)= VCG(1)#WO1) 
Where LM and VM are the longitudinal and vertical 


moments of the weight groups, respectively. 
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The designer is now given a screen display of phe 


total weight of lead, and is asked to input the desired 


amount of that total to be used as margin lead Cia 
This lead will arbitrarily be located at LOA/= 
longitudinally and D/=2 vertically. The remainder of the 
lead is considered to be trim lead (TL). The required 
LCG for trim lead is then calculated: 

LCGCTLI= CWCNSC)*LCBNSC -LM( Cond. A-1) -LM¢VL) 
mia wea 7 2 3/7 OL 

Similarly, the actual VCG(Lead): 

VCG(Lead)= (TL*¥6& +ML#¥R)/W( Lead) 

Where 6 feet above the keel is arbitrarily chosen 
for VCG(TL) to make placement in the circular cross 
section of the ballast tanks feasible. 

The designer is then asked whether the required 
LCG(SL) is feasible. If the response is affirmative, the 
program continues on to main ballast tank ae) 
Calculations. If not, the designer is asked for new LCG 
and VCG data for the weight groups and variable load. He 
is Prompted, in each case, with the previously entered 
LCG and VCG for each category. Of courses if reasonable 
adjustments in weight group locations cannot lead to 
feasible stability lead LCG (feasible being, generally, 
between the physical centers of the main ballast tanks), 
a re-appraisal of the pressure hull design and rough 
arrangements 15 necessary. 

When the lead has been satisfactorily located, 
calculation of required MBT LCG is accomplished 


(VCG(MBET) is arbitrarily fixed at R): 





aes 


LCG(NSC)= CLM(Cond. A-1) +LM(Lead) 
rain Vio a/7WCNSC ) 
LCG(MBT)= CWC Submerged)*LCBSUB 

~LCG(NSC)#WC(NSC) I/WCMBT ) 

The designer is again asked to decide whether the 
calculated required LCG is feasible, with Similar 
results in program flow depending on the aAnsueir. 
Feasibility of MBPT location is dependent on the desired 
geometry of physical locations for MBTs as envisioned by 
the naval architect. The continued graphic display of 
envelope and pressure full geometry is useful in aiding 
this decision process. 

The program now sums all vertical moments for 
weight groups, lead, variable load and main ballast to 
arrive at an overall VCG in the submerged condition 
(VCGSUB). An estimate of submerged stability may be 
calculated and displayed: 

BG= R-VCGSUP 

Note that the vertical center of buoyancy is 
assumed to be at R (D/2) for the submarine, which 1s 
basically a body of revolution. Any error due to 
buoyancy of appendages will be small and conservative, 
since off-center buoyancy items (Sail, eLew are 
generally above the geometric center of the envelope. 
Again, the designer is asked to accept or reject the 
calculated stability, and, again, rejection results in 


looping back to the input section of the program module. 


When all calculated values are deemed acceptable, 
the input data and results are written to an output file 
for record purposes, and the locations of each item’s 
center of gravity are plotted on the pressure full 
graphic display. The designer is now ready to proceed to 


mane calculation’ of the equilibrium palugon. 


igi 


THE EQUILIBRIUM POLYGON 


GENERAL DISCUSSION 


One of the most important sets of calculations in 
the submarine design process is the determination of the 
equilibrium polygon. Any submarine experiences changes 
in weight distributian during the course of operations, 
as items in the variable load are consumed. A system of 
tanks must be sized and located in Such a way that this 
change in weight distribution and amount may be 
compensated to allow the Submarine to remain in proper 
balance and trim. 

A detailed list of the variable load items to be 
considered may be found later in this chapter. In 
aggregate, the items compose major groups in the 
categories of provisions, ammunition, stores, fuel, and 
other consumable fluids. For a nuclear or other 
non-fossil-fuel submarine, the boat will always become 
lighter as items are consumed, and the campensation 


consists et veel ts) appropriate amounts) af variable 


ballast at the proper locations. In the case of a diesel 
or other fossil-fuel sttbmarine, a paradox exists for the 
designer. As fuel is consumed from the fuel ballast 
tanks external to the pressure hull, these tanks are 
compensated with seawater, since they are "soft" tanks, 
Tetiaiot es remain full of fluid to maintain a zero 
differential pressure across the tank structure. The 
seawater used for compensation has a higher density than 
the consumed fuel which it displaces. Consequently, the 
submarine actually becomes heavier as fuel is consumed. 
The tool which is used to calculate the required 
locations and capacities of the compensating tanks 15 
called the equilibrium polygon, which if a plot of 
weight versus moment (fore and aft) of the weight. 
Figure 7.1 15 an example of the equilibrium polygon far 
a submarine. The points within the polygon represent 
extreme operating conditions for the submarine in terms 
of the weight and moment required to balance the 
conditions. The solid lines forming the "polygon" are 
plots of the weight and moment added by progressive 
filling and emptying of the major compensating tank 
groups. These groups, which may be comprised of several 
tanks each (the aggregate weight and moment for the 
group is plotted), are the forward trim, auxiliary, and 
after trim tank groups. If all of the operating points 
for the submarine are contained within the polygon, the 


submarine may be safely compensated for all expected 
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SeopeeoneiecenGitions. [Tf not, the desiigner must change 
Poremboeattem(S), weight(s), or both, for the appropriate 
group(s), in order to expand the polugon to include all 
the conditions. 

The conditions to be calculated for polygon points 
are specified in the NAVSHIPS Technical Manual (NSTM) , 
wnaeter s27Y8 Cis i. The extreme conditions are more 
severe than an actual submarine would be expected to 
encounter, and represent a “worst case” situation for 
weight and moment balance. These extreme conditions may 
be thought of as representing the following operational 


situations: 


HEAVY #1: Following a short, fast patrol, during 
which no ordnance is expended, but all fuel oil is 


consumed. 


HEAVY #2: Same as the previous condition, except 
that only the fuel oil in the fuel ballast tanks 1s 


consumed. 


Lieem tieé LIGHT #2: Two variations following a 
stort patrol with no fuel consumed, but with all 


ordnance expended. 


Peeve QRWARD #ic¢ Following a patrol where only aft 


ordnance is expended, and all fuel oil is consumed 


from the forward fuel ballast tanks 


HEAVY FORWARD #2: Same as previous Geman Morin 
except that fuel from normal fuel oil tanks is also 


consumed. 


HEAVY AFT? Following a patrol with only forward 
ordnance expended, and with fuel oil from the aft 


fuel ballast and normal fuel oil tanks consumed. 


CONDITION Nt The "normal" full-load patrol 
GGnea2ttion at the beginning of a patrol. For 
submarines with FETs, the FRTs are specified as 


fully ballasted with seawater. 


CONDITION M: The same as condition Ns except that 
Submarines with FETs carru a full fuel load in the 
PEs (for sien slibmarines, this condition will be 


hoe teer. than Condition WN). 
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THE VARIARPLE LOAD 
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The items that comprise the Category referred to 
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the variable load are numerous and diverse ion Nacwre. 
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lireation of the items LS again aided bey binge 


guidelines and categorie: 


til 
fi 
G 
tT 
7 


-ified in tne MeaVSHIPS 
Technical Manual. 


Them variable load 1S Sttgqiyvided Into two 
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SUD—groups 5 fixed items, nec expected fae var 
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Soni ttecantiy during a patrol, and the ee Variable 
items, routinely expended as the submarine operates. The 
wt eG Partition of the variable load COnNSists fate the 


TOoOlloOWInGa items: 


CeewentbD EFFECTS: Actual weight of the submarine’s 
complement of personnel, clothing, and associated 


personal effects. 


SLEMs OR COMPENSATING WATER: The weight, for a 
fleet ballistic missile (FRM) submarine, of the 
submarine launched ballistic missiles (SLBMs) or 
the compensating water in their absence. SLEM 
compensating water exactly matches the weight and 
moment of the missiles themselves, so that this 


category may, indeed, be considered a fixed item. 


SANITARY TANKS: The weight of tanks and contents 


associated with waste disposal. 


OXYGEN CANDLES: The weight of these devices, used 
for the generation of oxygen under emergency 
conditions. These are generally a very small weight 


le ine 


PEeeotl IN SUMPS: That weight of lube oil 


Gemcalined in the oil sumeps of equipment. 


FIXED CLEAN FUEL OTL? That weight of oil carried in 
non-compensated tanks which are maintained 
Saori (1.2... Shield tankS on nttclear 


submarines). 


AS can be seen from Reference [iz], there are other 
items specified as fixed load, but they are minor in 
nature, and beyond the level of detail required for 
conceptual design. For example, the item "depth cantroal 
tanks" in the NSTM is not included. This omission will 
be explained in the discussion of fuel load subdivision 
in the following sectian. 

The items considered variable in nature are as 


follows: 


PROVISIONS AND STORES: The consumable foodstuffs 


aS hy 


(Provisions) and spare parts, paper Supplies, etc. 
(stores) that the submarine crew will use during a 


Patrol. 


REVITALIZATION OXYGEN: Oxygen used for life support 
(atmosphere control) during lengthy submerged 
periods. On newer submarines (particularly 
nuclear), a method of oxygen generation often 


obviates the need for oxygen banks. 


TORPEDOS, MISSILES, AND AMMUNITION: This category 
is self-explanatory. All expendable ordnance other 


than SLEMs is included. 


WRT TANKS: This item is the “water round torpedo" 
tank capacity, used to compensate for the presence 


or absence of torpedos in the tubes. 


RESERVE ELECTROLYTE: The electrolyte carried to 


replenish the submarine storage battery. 


FRESH WATER: For this thesis, this 1S a composite 


category, consisting of the total of potable, feed, 


and battery water. 


RESERVE LUPE OIL: Lube oil carried to replenish the 


Slime lube O11 during the patrol. 





FUEL?S Propulsion fuel for the submarine. See the 
discussion that follows for the sub-categories of 
fuel and their significance to the polygon 


computations. 


Fossil-fuel submarines present a variable load 
parodox for the designer. Fuel which is carried external 
to the pressure full (the majority of the fuel load “for 
most "conventional” submarines) must be replaced by 
seawater as it is burned, Since these external "fuel 
ballast tanks" (FBTs) are not constructed to withstand a 
high pressure differential (1.8., they are "SoT tT 
tanks"). The seawater which 1s used to compensate these 
tanks 15 heavier than the fuel it replaces. Thus, as the 
submarine consumes fuel, it gets heavier, rather than 
lighter. 

To provide a margin of safety at the base of the 
equilibrium polygon, some means mist be provided to 
remove weight as fuel is burned. In earlier designs; a 
"Safety tank" was provided, of sufficient capacity that 
when it was emptied, enough water weight was pumped 
overboard to compensate for the added seawater in the 
FRTs as fuel was consumed. The safety tank, then, was 
full at the beginning of the patrol, and would be empty, 
or nearly so; when all fuel had been removed from the 


FRTs. For purposes of this thesis, the safety tank 


concept was not used, since the weight and volume of the 
tank and it’s contents is a non-productive load item. 
However, a means to compensate for the FPT weight 
addition 185 necessary, 

The method chosen for implementation in this thesis 
is the provision for a variable fuel oil (VFO) tank. The 
VFO tank is internal to the pressure full, and does not 
require seawater compensation as fuel is consumed. fn 
actual practice, this tank may be ballasted and treated 
as a dual-purpose tank, acting as an additional trim or 
auxiliary tank. For computational simplicity within the 
Program module, however, the VFO tank is treated as an 
uncompensated clean fuel o11 tank. An operational mode 
of employment is assumed whereby the fuel in the VFO 
tank 1S consumed first during the patrol, before the 
FPTs are utilized. Thus, the size of the VFO tank must, 
at a minimums, be sufficient to compensate for the 
difference in fuel and seawater densities for the FBT 
Capacity. This dictates a minimum capacity of 23% of the 
total fuel load. The designer is asked to input the 
percentage of tatal fuel to be allocated to the VFO tank 
(or tank group, as arrangements may dictate). I f 
sufficient internal volume is available, it may be 
advantageous to use a figure greater than 23%, in order 
to provide an adequate margin at the base of the 


Polygon. 





CALCULATING THE POLYGON 


The program module execution begins by reading the 
weight group values, submerged LCB, and LCG of the main 
ballast tanks from data files generated by the weight 
estimating and balance modules. The total Ilnmad to 
submerge the submarine is calculated: 

Load to Submerge= Submerged Displacement- Condition 

A Weights. 

The total variable load weight is displayed, and 
the designer is prompted, item by item, for the weight 
and LCG of each variable load item. As each weight = and 
LCG are entered, the balance of the weight in the 
variable load account is displayed. Items which are 
commoniy distributed in two sub-items, fore and aft, are 
input as two separate weights and LCGs. Examples of such 
items are torpedos, missiles; and fuel. The designer 
should have given some thought prior to this input 
sequence to the subdivision of the variable load. The 
envelope and pressure full are displayed at the top of 
the screen during this stage of the module, to 
facilitate estimating the LCG for each item. 

The last item to be entered is fuel. The balance of 
the variable load account at this point is assumed to be 
all fuel (obviously; the balance should be zero for 4a 


non~fossil-fuel submarine). nice aot acl fuel load is 





displayed, and the designer is asked to enter the 
fraction to be allocated to the VFO tank. He is then 
asked to provide LCGs for the forward and aft FBTs and 
the VFO tank. 

With the input of all items completed, the program 
displays a variable load summary of all item weights and 
ECGs’. If the summary is satisfactory, the program 
Proceeds with calculations; ie Went the designer is 
returned to the input sequence, where the prompt for 
each item now includes a display of its current value. 
If this value is to remain unchanged, the operator 
Simply hits the "RETURN" keys; otherwise, a new value may 
be specified. This iterative input procedure is repeated 
until the designer is satisfied with the variable load 
Summary. 

When the aiunput of all variable load items is 
complete and accepted by the operator, the program 
calculates the water to balance, and the associated 
moment of that water, for each of the nine specified 
POlygon points. Note that conditions M and N will be 
identical for nuclear submarines. Condition M is only 
defined for submarines with FPTs. 


A computation of one of the equilibrium conditions 


would Proceed in the following manner. Prior to 
computing the individual conditions, the arms and 
moments for each of the variable load items are 


determined: 


Item Arm= Submerged LCR —Item LCG 

Item Moment= Item Weight *Item Arm 

In addition, factors of Proportionality are 
determined for light and heavy seawater. For each 
condition, the sum of all variable load items, in the 
fractional amount specified by the NAVSHIPS Technical 
Manuals 15 computed. Similarly, the total moment of the 
specified items is determined. The required water to 
balance and moment to balance may now be calculated: 

Water to Palance= Load to Submerge -MPT -Sum of 

Item Weights 

where MET= Main Ballast Weight 

Moment to Balance= -MBT Moment -Sum of Item Moments 

If required, the water and moment to balance are 
ieee lied by the proportionality factor for Ddight or 
heavy seawater. Each of the nine Conaderans 1s 
calculated in a Similar fashion, resulting in arrays of 
weight values in tons and moments in foot-tons (with 
negative moment values indicating an aft moment). 

Several "housekeeping" chores are performed at this 
Point in the module. The display of envelope and 
pressure hull, and associated data base attributes, are 
filed under the part name MIT.HULLOUT, for later recall 
to plot desired views of this graphic screen. The 
maximum values of water and moment to balance are 
examined to determine the required scale for the polygon 


Plot. A new part is activated under the name MIT.POLYOUT 





for graphically displaying the polygon. Depending on the 
scale required, one of several pre-constructed "form" 
drawings 1S activated, providing the AXES and 
appropriate scale information for the polygon under 
consideration. 

The program now inserts and labels each of the 
Perm rirecOnad2t1OoOn points on the graphic display. 
This gives the designer a clear display of the limits 
required for his Polygon design. An important 
consideration at this point is whether an adequate 
margin remains at the base of the polygon. Generally, if 


the lowest weight value 1s less than 15 tons from the 


abscissa, the designer should seriously consider 
re-structuring the Variable loads as necessary, to 
achieve this margin. This may be accomplished by 


returning to the input section of the module. An example 
of such a change would be to place more of the fuel 
load, volume permitting; in the VFO tank, for a 
fossil-fuel submarine. 

The operator is now asked to provide the LCGs of 
the Forward Trim, Auxiliary, and After Trim tank groups. 
He is then prompted for the capacities; in tons; Or 
each of the tank groups. The program determines the arms 
and moments of each group in the same fashion previously 
delineated for the variable load items. 

Lines are now inserted on the display to represent 


the successive filling of Forward Trim, Auxiliary, and 





After Trim tanks, followed by the successive draining of 
the tanks in the same order. These lines form the 
Polygon, and all equilibrium points must fit within its 
boundaries if the tank locations and capacities are 
adequate. Obviously, a POlygon whose boundaries 
substantially exceed the extent of the equilibrium 
Points, while safe, represents a waste of internal 
volume, and should not be accepted. 

The designer now has the option to accept or reject 
the polygon, as Plotted. If it 15 rejected, the polygon 
boundaries are erased, and the designer may specify new 
tank group capacities, LCGs, or both, after which a new 
Polygon is plotted. This process may be repeated as many 
times as are required to achieve an acceptable polygon. 

When the polygon results are satisfactory, the 
pPolugon is plotted on the electrostatic plotter, the 
part MIT.POLYOUT ais deleted, and the variable load 
summary, tank locations, and tank capacities are written 
to an output file for a permanent record of the final 
input data set. One iteration of the submarine 
conceptual design process, as supported by this thesis, 


15 complete. 


SUMMARY AND 


RECOMMENDATIONS FOR FURTHER RESEARCH 


SUMMARY 


This thesis constitutes a first exploration of the 
use of interactive graphics for conceptual design of 
submarines. The resulting product is useful, flexible, 
and accurate for conceptual design activity. It is by no 
means a finished product in the commercial sense, but 
Provides a tool to build on and to use confidently for 
such applications as student design projects. 

The power of real-time graphic display of the 
design during its creation 15 a major factor in the 
worth and potential of the program. Many hours of 
laborious mantial dvat ting are avoided by the 
high-resolution graphic output. In addition, a user with 
even a rudimentary knowledge of the Computervisioan 
system can go far beyond the capabilities of the basic 
software package, thus greatly enhancing its value as a 


design tool. 





View Progream 15 highly interactive; the Liser 
controls all important design decisions. The designer 
must have a solid knowledge af naval architecture, which 
1S appropriate for any ship design program. Subjective 
design ygudgements must remain in the province af the 
user for a credible result, and that philosophy has been 


rigorously adhered to. 


RECOMMENDATIONS FOR FURTHER RESEARCH 


There are nearly unlimited opportunities For 
continued work on this project, and ain this general 
area. This thesis has scratched the surface in an area 
ripe for research and development. 

The first obvious areas encompass those portions af 
the design process which were excluded from the current 
project. These include structures, propeller design and 
detailed interior arrangements. 

Submarine structural design could benefit greatly 
from treatment in the CAD environment. Visualization of 
complex structural relationships is intrinsically 
Clearer with the use of real-time graphic display. In 
addition, the use of powerful finite element packages, 
with automatic mesh generation, 1s possible on many CAD 
systems, IWewiding Computervision. These Packages 
display stress levels in color, and will magnify 


displacements on command for rapid visual analysis. 


Arrangement of interior details is feasible on 
these systems, including automatic interference checking 
eae lavient of Piping and electrical Gistri1 bution 
networks. 

Propeller design may be facilitated by graphic 
display of the flow regime during analysis. Much 
analytical work remains in this area, and research into 
links with main-frame systems and large data hases is 
also needed. 

In addition, the translation of this package into 
the Fortran-S language, with attendant PYroOodgranming Of 
linking and loading files, would greatly improve 
execution speed and I/O flexibility. The creation of 
customized graphics commands would add more power” and 
enhance the ergonomics of the package. 

Further work in the subject areas of weight 
estimation, hydrodynamics, and animated display Of 
control system performance is required to derive the 
maximum benefit from submarine CAD. All = of these 
subjects have been separately addressed in other 
software Packages, but would be valuable additions to an 
integrated design program. 

With increased interest and design volume in the 
submarine design field, research funding 1s bot ti 
necessary and clearly justified by the potential 
improvements to the speed and accuracy of the design 


Process. The overall area of computer aided design Or 





submarines can support major research projects in the 
future, benefitting both academic and commercial design 


metivities. 
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MIT.&BCD.WEIGHT 
3-28-84 13337244 FUTIL 6.18 


I!WEIGHT ESTIMATE FOR 13.461 DESIGN PROJECT: 
eo 


3!TYPE: ATTACK 

4!PROPULSION:s MODIFIED WANKEL ELECTRIC 
S!BATTERY TYPE NICKEL CADMIUM 
6!SPRINT KW-Ht 3300 

T! ENDURANCE KW=Hs: 8000 

BS!BATTERY WEIGHT: 165.631 

OS! ENGINE HORSEPOWER: 4500 

JO!ENGINE & ASSOC. EQUIP. WEIGHT: 22.86 
Pt!MOTOR & ASSOC. EQUIP. WEIGHT: 32.25 
l2!DISCRETE ADDITION TQ GROUP 2 WEIGHT: 80 
A3IMAX. OPERATING DEPTH IN FEET 700 
14! 

TIS5!LEAD FRACTION: Q.1 

I6O!VARIABLE LOAD FRACTION: 0.089 
17!RESERVE BUOYANCY FRACTION: 0.125 
18!FREE FLOOD FRACTION: 0.08 


iS! 

Zul! 

21 !CATEGORY WEIGHT 
22! 

23!GROUP 1 493.318 
24!GROUP 2 300.741 
25!GROUP 3 17.4418 
26!GROUP 4 95 
27!GROUP 5 137.043 
28!GROUP 6 lites 
29!GROUP 7 Se 
30!COND. Al 1245.84 
Sit LEAD 124.584 
32!COND. A 1370.43 
SS!) VAR. LOAD 133.884 
34!NSC 1504.31 
35° MST | 88.039 
36!SUB. DISP. 1692.35 
BeerRee FLD. 135.388 
38!ENV. DISP. (827.74 





MIT. &BCD.GEQMOUT 
3-28-84 13238256 FUTIL 6.18 


b!GEOMETRY OUTPUT FOR 13.46! DESIGN PROJECT: 
a. 


oe 

4'tLENGTH OVERALL’ 139.947 

S!MID=BODY LENGTH? 74.9465 

6!OVERALL PRISMATIC: 0.798368 

7! LENGTH/DIAMETER? 5.1832 

8!SUBMERGEN LCB (AFT OF FP)! 66.2489 
Q!DRAFT @ NORMAL SURF. COND.? 21.9585 
1O!LCB @ NORMAL. SURF. COND. ! 66.3329 
rele 

}2!RECAP OF INPUT DATA? 

3)! 

I4tENVELOPE DISPLACEMENT #1827. 74 
IS!tENTRANCE? 25 

16!RUN : 40 

UTIDTAMETER? 27 

IS!FWD PRISMATIC: 0.6 

IS9'tAFT PRISMATIC: 0.55 





MIT. &BCD. OFFOUT 
3-28-84 131%38:56 FUTIL 6.18 


1 !LENGTH WRT FP OFFSET 
2!0 0 

JUS 1.79298 
411.5 3.41525 
5!2.99946 5.08919 
6!4.99839 6.77129 
7!6.99732 8.1157 
8!13.9946 11.36/74 
9120.92 13.1284 
10!27.9893 ; 13.5 
11!34.9866 (325 
12!41.984 13.5 
13!48.98)3 32> 
14155.9786 (225 
15!62.9759 i335 
16!69.9733 13.5 
17! 76.9706 13.5 
18!83.9679 13.5 
19$90.9652 13.5 
20!97.9625 13.5 
21!104.96 13.334 
22 aie. 957 12.444 
23S. 255 10.7078 
24!1125.952 8.07667 
25!132.949 4.51632 


26!139.947 0 


Sy 


MIT.&BCD. POWER 


4= 9-84 10:13:59 FUTIL 6.18 


1'SPEED AND POWER RESULTS FOR 13.461 DESIGN PROJECT! 
2! 
3'NON=PROPULSION LOADS (KW)? 135 


4!PROPULSIVE COEFFICIENT s O.77 
( 

ri 

T!SPEED (KT) SHP TOTAL KW 

8! 

9! | 0.615861 135.459 
1O Nee 4.61783 138.444 
Pits 15.0496 146.222 
12! 4 34.8408 160.981 
ist) 5 66.8572 184.855 
14! 6 113.92. 219.95 
15! 7 178.818 268.345 
16! 8 264.312 332.097 
17! 9 373.141 413.251 
18! 10 508.023 513.833 
9! 1] 671.663 635.859 
DO Lee 866.75 78) .335 
2i!i 13 1095 .96 952.256 
22! 14 1361.95 1150.61 
23! 15 1667.38 1378.37 
24! 16 2014.9 1637.51 
25! 17 2407.13 1930 
26! 18 2846.71 2257.79 
a7! 19 3336.24 2622.83 
28! 20 3878.34 3027.08 
29! 21 4475.63 3472.47 
30! 22 5130.68 3960.95 
31! 23 5846.11 A494. 44 
32! 24 6624.47 5074.87 
33! 25 7468.37 5704.17 
34! 26 8380.38 6384.25 
35! 27 9363.07 7117.04 
36! 28 10419 790 4. 45 
37! 29 11550.7 8748.37 
38! 30 12760.8 9650. 74 
39! 31 14051.8 10613.4 
40! 32 15426.3 11638.4 
41! 33 16886.8 12727.5 
42! 34 18435.8 13882.6 
43! 35 20075.9 15105.6 
44! 36 21809 .6 16398.4 
45! 37 23639.5 17762.9 
46! 38 25568 19201.1 
47! 39 27597.7 20714.6 
48! 40 29 731.2 22305.5 





MIT. &BCD. BALOUT 


5- 


6-84 


[6254237 FUTIL 6.21 


7h = 


}!RESULTS FROM BALANCE MODULE 


2! 
3!GROUP WEIGHT LCG VCG 

4! 

5! 360.354 110 15.5 

6! 2 623.431 120 1U 

(3 00,3814 I12 15 

Bi 4 105 60 20.5 

9! 5 243.99 110 9.5 
10! 6 195.5 10 22.5 
(172.90 90 10 

12! Al 2218.09 106.19 13.902 

13! LEAD 221.809 16.889 8.1414 

14! A 2439.9 103.53 13.378 

15! V LD 183.643 9% 16 

16! MBT 327.943 101.39 15.5 

17! 

I8!MARGIN LEAD (TONS)! 50 @VCG=N/2, LCG=LO0A/2 
IO!STABILITY LEA) (TONS)? 171.809 AVCG= 6 FT 
20!STABILITY LEAD REQUIRED LCG: 67.1016 
21! 
22!NSC (TONS)! 2623.55 ¥LCBt 103.145 
23! 
24!SUBM (TONS)? 2951.49 @LCBt 102.951 
25! 
20!SUBMERGED STABILITY (BG)? 1.72272 FT 


5= 7-84 19:24:02 FUTIL 6.21 


1!PQLYGON SUMMARY OUTPUT FILE 
2! 


3! 

4! ITEM 

5 

O'CREW AND EFFECTS 
T!SLBM7S OR COMP. WATER 
BISANITARY TANKS 

SO EUBE OLESEN SUMPS 
IO'FIXED CLEAN FUEL OIL 
1itw*e*nTQOTAL FIXE!) ITEMS **x% 
12!PROVISIONS AND STORES 
IS!REVITALIZATION 02 
14¢TORPEDOS IN FWD ROOM 
1S'tTORPEDOS IN AFT ROOM 
I6!TACTICAL MISSILES FwWN 
17!TACTICAL MISSILES AFT 
ISIwRT TANKS 
IS!RESERVE ELECTROLYTE 
2O! TOTAL FRESH WATER 
21!RESERVE LUBE QIL 
P2'FUEL BALLAST FWD 
23{FUEL BALLAST AFT 
2a VAR VADEE FUEL OLE 
25! 
26! 


ECG 


55 

O 

G2 
90 

0 
tees 3 
1S 
88 
45 

O 

O 
116 
30 
58 
100 
95 
26.3 
Lies 
45 


27!VARIABLE BALLAST TANK DATA? 


23! TANK GROUP 
29!'FORWARD TRIM 
30! AUXILIARY 

BIG ARIER TR1M 


ECG 
oD 
1 
105 


WE IGHT 


CAPACITY (TONS) 
25 
100 
55 





Almeenorx. it: 


SevVRcE CODE LISTINGS 


Bea 


MIT.&8CD.DESIGN 
5- 7-84 13250200 FUTIL 6.21 


fore WC17).PC17) . FCS) RS 21) CIRC 21) .LX(8) .RXC3) .XPO0S(30) .SAC21) <# 
Otel 5S(21).HPC41). KW(41) UC21) .X029).Y¥( 29) 

3!READ CENTER DESIGN NAME? ) A&NAME 

4!#WEIGHT CONTINUE 


5!N=| 

6!XEPEAT 

i} N=N+/ 

3! IF (N.EQ.3.0KR.N.2E0.5) GOTO LOOP 

9! PRNT DO YOU WANT TO SPECIFY (S) OR CALCULATE (C) GROUP (N] WEIGHTS? 
10! READ AR 
li: IF (&R.FEO."S") FCN) =I 
i2! IF (&@RCEO."C") FCN) =0 
13! IF (FC(N).EQ.0) GOTO LOOP 
14! PRNT ENTER GROUP [EN] WEIGHT (TONS). 
15! READ WT 
16! W(N) =NT 


17!#LOOP CONTINUE 

ISO!UNTILON. EQ. 7) 

IStREAD( MAX. OPERATING NEPTH IN FEET: )ZMAX 
20!P(1)=.215541. 606225E-4 #ZMA K 

21 !QEADC RESERVE BUOYANCY FRACTION? )P(13) 
P2'READCLEAD FRACTION: )P(9) 
O3'!READ(VARIABLE LOAD FRACTION: )PCTI) 
24!READXC FREE FLOON FRACTION: )P(15) 

25!PRNT NUCLEAR (1). DIESEL (2). OR WANKEL (3) PROPULSION? 
Z26!READ PM 

27!PRNT ATTACK (1) OR FRM (2)? 

23!QEAD T 

2O!TF (F(2).EQ.1) GOTO LABI 

30!GOSUB GROUP2 

31!#LAB! IF (F(6).EQ.1) GoTO Al 

32!GOSUB GROUP6 


33!GOTO Al 

34! #GROUP2? 

a5! READ (SHAFT HORSEPOWER: )SHP 

30! IF (PM.NE.1) GOTO BATT 

37! K=.O03*SHP/100041.13 

33! WI=(K* LES *SHP) /(LGCSHP )) **5 

39! A(2)SWI/ 2240 

40!RTNSUB 

41!#BaTI CONTINUE 

42! PRNT BATTERY TYPE: NICAN (1). VARTA (2). OR TRIDENT (3)? 
43! REAL) BT 

44! REANC TOTAL ENGINE HORSEPOWER: )ESP 

45! IF (PM.E9.2) KE=. 0096 

46! IF (PM.EQ.3) KE=.00508 

a7! WE=KE*ESP 

4c! WM=,00645*SHP 

49! READ(TOTAL KW-H STORAGE FOR SPRINT: )KSP 
S50! READ(CTOTAL KW-H STORAGE FOR ENDURANCE RATE? )KEND 
51! Dewees) C5=35.5610 

52! Le eee. 2) CS=22. 353 

53! le 4 B1.6E9.3) CS=12.1925 

54! B=KSP/CS 

B5! Mee oieeo. I) CE=46. 3 


DO. fepetnhe eo. 7) CE=47.074 





BO 


pi! Peer EO 47) CE=25.401 2 
53! Bl1=KEND/CE 

59! Pe C81 .Gl oR) B=B 1 

60! WT =WE+WM+B 


ol! PRNT ENTER DISCRETE ADJUSTMENT TQ GROUP 2 (0 OR # OF TONS). 
62! READ DA 

63! W(2)SWT4+DA 

O4!RTNSUB 

65! #GROUP6 

66! HEAD(NUMBER OF PERSONNEL: ) NP 

6/i! W(6J=2.3*NP 

63 !RTNSUB 

69! #A] 

TO!IF (PM.LEQ.1) P(3)=.045 

TILTF (PMCEQ.2.0OR.PM.EQ.3) P(3)=.014 

7292 (4)=.06 

Wee le {(T.EO.!) P(5) =. 11 

f4'!1F (T.EQ.2) P(5)=.095 

PoliF (7.EC.1.AND.PM.EO.1) P(/)=.03 

f6O'IF (T.EQ.1.AND.PM.EQ.2) P(/)=.04 

V7T!ITF (T.EQ.!1.AND.PM.EQ.3) P(7)=.05 

fMSEIF (T.FO.2) P(/J)=.14 

TON FW=(P(1)4P01)*P(9))/4C1-PC1I1))+P03)4P(5) 
Solte. (FP (4) .EQ.0.ANOD.F(7/).EOQ.0) GOTO CASE! 
SI!iF (F(4).EQ.1-AND.F(7).E0.9) GOTO CASE2 
S2!IF (F(4).EQ.0. aND.F(7).FO.1) GOTO CASE3 
O3!1F (F(4).EQ.1.AND.F(7).EQ.1) GOTO CASE4 
34! #CASE] 

55! AlL=(N(2)4W(6))/ 01] -FW=P(4 =P /)) 

36! W(4)=P (4) %Al 

af! W(/)=P(/) «AI 

83!GOTO TOTAL 

sO! e#CASE2Q 

90! Al=(N(2)4N(4)4N(6))/S (1) 7 FNMP(7)) 

Q1! W(7)=P(7)#Al 

92!GOTO TOTAL 

93! ACASE3 

94! Al=(W(2)4N(6)4N(0/))7(] —FW-P(4)) 

95! W(4)=P(4) A] 

96!GOTO TOTAL 

O7!#CASE4 

9d! Al=(W(2)4N(4)4N(60)4W0/))701 -FW) 
99! # TOTAL 
100 !W(3)=P(3)*#Al 

IOLIW(5)=P(5) #A] 
IO2'W(d) SA] 
1IOZ!IW(9)=P(9) ¥A] 

LO4tw( 10) =A] +W(9) 
LOSIH (TOA, HSRC ITO) C1 =P C11 )) 
IOSIWC 11) =PC11 )*#W0 12) 
IOTIN(]) =PC 1) W112) 
JOGIWw(13)=P013) "WwW (1 2) 
IOQ TACT 4) =AN(12)4W(13) 

LIVIN (15) =P 015) #N (14) 

LITIWNC 16) EN015)4W014) 

Piet) CALCUMATIONS COMPLETE... 

113! PANT 

114!PRNT 

bi5¢<* OUTPUT SECTION 

Diente.) SPata Ty acy 





a0 — 


(T.EQ.2) &T="BALLISTIC MISSILE* 
(PMZEQ. 1) &PM="NUCLEAR# 

(PM.EQ.2) &PM="DIESEL ELECTRIC" 
(PM.EQ.3) &PM="WANKEL ELECTRIC" 
(DA.NE.O) &M="MODIFIED * 

(DA.EQ.0) &M="4 

(BT.EQ.1) &BT="NICKEL CADMIUM" 
124!1F (BT.EQ.2) &Bl="LEAD ACID (VARTA) 
I25'IF (BT.EQ.3) &BT="LEAD ACIO (TRIDENT)" 
126!PRNT WEIGHT ESTIMATE FOR [(&NAME]J? 

127! PRNT 

123!PRNI TYPE: (AT) 

129!PRNT PROPULSION? (8&M)(&PM) 

I30!IF (PM.EQ.1) GOTO SHAFT 

131!PRNT BATTERY TYPE: (&BT] 

132! PRNT SPRINT KWH? (KSP] 

133!PRNT ENDURANCE KWeH: [(KEND]) 

134!PRNT BATTERY WEIGHT: [(B]) 
135!PRNT ENGINE HORSEPOWER? 
136!PRNT ENGINE WEIGHT: [WE] 
I37!PRNT MOTOR WELGHT: (WM) 
PRNT DISCRETE ADDITION TO GROUP 2: 


LI7tIF 
es! LF 
L1IOITF 
Poo! 1F 
PZ131F 
22! 1F 
Pes IF 


(ESP) 


(DA] 


13u!IF (DA.NE.O) 


I39!#SHAFT CONTINUE 


140!PRNT 
141!PRNT 


l42:TF (FC6).EQ.1) PRNT NUMBER OF PERSONNEL? 


143!PRNT 
144!'!PRNT 
145! PRNT 
146! PRNT 
147! PRNT 
l4ad! READ 
149!PRNT 
150! PRNT 
151 !PRNT 
152! T=! 


{SHP ] 
(FT): 


SHAFT HORSEPOWER: 
MAX. OPERATING DEPTH (ZMAX]) 

(NP J 
LEAI) FRACTION: [(P(9)] 

VARIABLE LOAD FRACTION: (PC 11)] 

RESERVE BUOYANCY FRACTION: (P(13)] 

FREE FLOOD FRACTION: [(P(15)] 

Ht “<RETURN> TO!) CONT I NUE. 

ACONT 


CATEGORY WEIGHT 


ia REPEAT 


154? PRNT 


155! l=I +1 


GROUP [(T} CWC) ) 


I56!UNTIL (1.EQ.8) 


15 /!PRNT 
158! PRNT 
159! PRNT 
160! PRNT 
lol !PRNT 
152!PRNT 
163!PRNT 
164!PRNT 
165!PRNT 
166! PRNT 
167!PRNT 
1S65!REAN 


(W(3) J 

(W(9)] 

(WC10)] 
(WC 11)] 
(W(12)) 
(W(13)]) 
(WC14)) 
(W(15)] 
(W(10)] 


COND. Al 
LEAN 

COND. A 
VAR. LOAD 
NSC 

MBT 

SUB. (WLSP. 
FREE FLD. 
ENV. ODISP. 
(R)? 


SATISFACTORY (5) OR RECALCULATE 


&NECIN 


ieee AWECID. EO. "Rh" IeGOTO WEIGHT 


170! PRNT 


172 PRNT 


PLEASE WAIT WHILE OUTPUT 


COMPLETED TRE ENTIRE DESIGN PROGRAM, "MIT. WEIGHT" 


AND PASS FILES ARE WRITTEN. 
Pee ete ew ElGd th CALCULATIONS WILL BE. IN FILE “MIT.WEIGHT". 


THE OUTPUT FOR 
WHEN YOU HAVE 
AND THE OUTPUT 


1/3!PRNT FILES FOR SJBSEQUENT MODULES MAY BE OBTAINED IN HARN COPY BY USING 


1/74!PRNT 


THE COMMANI) “PRINTLP"FILFENAME!'%, 


Pio OPENS.’ Vv IT Ai bOnT" 


a. ke 


ree 


| 77!'&S‘{P=SHP 

l/s! 8ESP=ESP 

1}79!'!&KSP=KSP 

130 !&8KENND=KEND 

IS 1!&ZMAX=ZMAX 

182 !&NP=NP 

183!&8BB=B 

184!&8WM=WM 

1 35!8WE=WE 

1S36!'&DA=DA 

ly7!&RLF=Pp(9) 

|} 388!aAVLF=PC 11) 

159 !8RBF=P(13) 

IS9O!&SFFF=P (15) 

IS T'AWIl=an( 1) 

192! 8W2=W (2) 

193!'&8W3=W(3) 

194 !8W4=nw(4) 

195! &8wW5=W (5) 

IS6!'!R&W6=N (6) 

1O7!8W7=W( 7) 

198! 8WN8=n (3) 

199 'tQAWwO=sn(9) 

200! 8W10=N( 10) 

QO1'S8wil=an (il) 

20218WI2=N( 12) 

2U3'8W13=N(13) 

204! 8WI1 4=N (14) 

205! 8W15=NW(15) 

2UO'BWISO=N(16) 

2U 7! 8X="WEIGHT ESTIMATE FOR “+8NAME4" 38" 
203 !GOSUB QO 

2U9!WRITEF 2.8P 

21IOF8X="TYPEs "+8T 

211!GOSUB Q 

212!8X="PROPULSION? "+&M+&PM 

ms ie) (PM. EQ. 1) GOTOUNUKE 
214!GOSU8 0 

Plo Ge="OATITERY TYPE: e@+a8T 
Z216!GOSUB Q 

21 7tAX="SPRINT KWH: "+&8&KSP 
ZI13!GOSUB 0 

Z2IDtAX="ENDURANCE KWeH: "+&KEND 
220'GOSUB 

221 !AX="BATTERY WEIGHT: "+&BB 
222!GOSUB Q 

AE23'aAX="ENGINE HORSEPOWER! “+8ESP 
224!GOQSUB 

225!8X=S"ENGINE & ASSOC. EQUIP. WEIGHT: "+8WE 
226!GOSUB 0 

POT!RAXS"MOTOR A&A ASSOC. EQUIP. WEIGHT: “+8WM 
P26!GOSUB 

PAQGITF (NA.EN.J) GOTO NUKE 
A230'AXS"NTSCRETE ADDITION TO GROUP 2 WEIGHT: "4&DA 
A31!I!GOSUB O 

232'e#NUKE CONTINUE 

233!AXS"MAX, CUPERATING NEPTH IN FEET: "“#&ZMAX 
234!G50SUB 0 

P35!WRITEF 2,&P 

PSO*RV="LL EAI) FRACTION: "+44LEF 





237'GOSUB O 


238!8&X="VARIABLE LOAD FRACTION? 


239!GOSUB 0 


Oe. 


UtAN EF 


240! &X="RESERVE BUOYANCY FRACTION:S "+8&RBFEF 


241 !GOSUB 0 


242'&X="FREE FLOOD FRACTIONS 


243!GOSUB 0 
244!WRITEF 2,8P 
245!WRITEF 2,.8P 


240! 8X="CATEGORY WEIGH ** 


24 /'GOSUB O 
243!WRITEF 2,&P 


249!8X="GROUP 1 US | 
250!GOSUB Q 
251!8X="GROUP 2 "+8W2 
252!GOSUB 0 
253'&X=s"GROUP 3 "+3W 3 
254!GOSUB 0 
255!8X="GROUP 4 "+ &EW4 
256!GOSUB 0 
257!8X="GROUP 5 + RWS 
253!GOSUB 0 
259 '8X="GROUP 6 4+ BW6 
260!GOSUB 0 
261!8&X="GROUP 7 + KW 7 
262!GOSUB 0 
263!8X="COND. Al "+8Wo 
264!GOSUB 0 
265! &8X="LEAD "+&8W9 


266!GOSUB 0 


26 7!&X="COND. A "+8W10 


268 !GOSUB 0 
269! &X="VAR. LOAD ora || 
270!GOSUB 0 


271'8X="NSC "+8Wwl2 


272°SOSUB () 


2738&X="MBT "+ BWt3 


2/4!GOSUB 0 


2foe6r2—"SUB. DISP. "+&Wla4 


2/6!GOSUB 0 


Bo OAS" REE FLO). "+ &Wwi5 


2/3!GOSUB 0 


279'8X=S"ENV. DISP. "“+8W16 


2380!GOSUR 0 

281!GOTO DONE 

23250 CONTINUE 
2o3!WRITEF 2,&X 
234!RTNSUB 

285! #DONE CONTINUE 

230 °OPENW 3,"MIT.PASS1" 
OSlIWRITEF 3,awl2 
2ZESsrwRITEF 3.6015 
2oU!WRITEF 3.8&Wlo 

290 !WRITEF 3.4NAME 

POLIT OPENAN 4,¢MIT. BALPASS" 
2U2!WRITEF 4, AW] 
PQO3!NRITEF 4, AW? 
29O4!WRITEF 4, &W3 
2OS!IWRITEF 4, &AW4 
JOASIWRITEF 4, 4A45 


"48 FFF 





2O7T!IWRITEF 4, &W6 
293!WRITEF 4. &W/ 
299!WRITEF 4, &W3 
300!WRITEF 4, &W9 
301!WRITEF 4, &WI10 
302!WRITEF 4, &Wtl 
303!WRITEF 4, AWl2 
304!WRITEF 4, &wWt3 
305!WRITEF 4, &W14 
306!WRITEF 4, &Wi5 
307!WRITEF 4, &WI6 


303!0PENW 5, MIT. PASSPH2" 
309!WRITEF 5,8W10 ~ 
310!WRITEF 5,8NAME 
31T!RUN NEW MIT.G3 
312!END 





— 84 — 


MIT.&BCD.G3 
fee O4 IN t08226 FULIL 6.23 


Pomme otel).ClRAC(21) .LX(8) ,RX(8),XP0S(30).SA(21) .LS( 21) <# 
2'>,HP(4)),KW(41) ,—DUC 21) .X(132),Y(132) ,SX(21) 
3! 4#GEQOM 

4'QPENR 1.°MIT.PASSI" 

5!READF 1,4W 

6'!DNSC=8W 

T!READF 1,8&W 

8!DFF=aw 

S!READF 1,&W 

1O!DENV=&wW 

tT !READF 1,&NAME 

}2!PI=3. 14159265 

13! #BRANCH CONTINUE 

14'READ(ENTER LENGTH (OF ENTRANCE? )LF 
IS!READCENTER LENGTH OF RUN: )LA 
T6'READCENTER HULL DIAMETER? )D 

17!READCENTER FOREBODY PRISMATIC COEFF.: )CPF 
I8!READ(C ENTER AFTERBODY PRISMATIC COEFF.: )CPA 
IS! L=LFeLa 
20!LL=LF 

21! INTV=L/20 

22!R=D/2 
23!NFH71.04 77 «CPF #*5465.8107*CPF 224-483. 3037*CPF 42%*34587.413 / «CPF ex2<# 
24!>—281. 7224*CPF+49 5876 
25! NA=379.65468CPA ##5-938.4 708 *CPA &* 44944, 8853eCPA &*3-471 ,.0872eCPAex2<# 
26!>4119.1465*CPA-11.3454 
27!LS(1t)=0 

28!LS(21)=L 

29 !J=s2 
30!REPEAT 
31!LS(J)=(CJ-1)*INTV 
B2Uxe=LS (J) 
33!GOSUB RAD 
34!RS(J)=RANS 
35!SA(J)=PI*RS( J) **2 
36!DNUC I) =SA(J) 
37! JjJ=J+] 
38 !UNTIL(J.EQ.21) 

39!GOSUB INTG 

40! VENV=35*DENV 

4) !NV=VENV=$ANS 

42'LMB=DV/(PI *«R x2) 

43!LL=LF+LM3 

44! LOA=LF+¢LMB+La 

45! INTV=L0A/20 

46!LS (21) =L0A 

47!1LS(1)=0 

48!J=2 

49!REPEAT 

5OILS (J) =(J-1) «INTV 

51! XLaels(J) 

52!GOSUB RAD 

53!RS(J)=RANS 

S54!CIRC( J) =2*PI*kRS( J) 

5S5!SACJI=PI*®RS (J) #2 

56! DUC IV=CIRC(S) 





Sen 


57! Jase) 

S3!UNTIL(J.EQ.21) 

59!Gosu8 INTG 

O0!SF=ANS 

61! VCYL=PI*LOA®Re¥2 
62!'CP=VENV/VCYL 

63!J=2 

64!REPEAT 

65!DU(C J) SLS (J) *SACJ) 

66!J=J+] 

6O7!tUNTIL(J.EQ.21) 

68!GOSUB INTG 

69 !LCBSUB=ANS/VENV 

70!DR=.5*R 

71 !DN=DNSC+.8*DFF 

72! LO=DN-.005*DN 
73'HI=DN+.005*DN 

14!PRNT DNt (DN) 

TSIPRNT LO: (LO) 

7T6!PRNT HI: (HI) 

7T7!#LCBLP CONTINUE 

T3!PRNT ITERATING ON NSC DRAFT @ T= [R*#DR]) 
719! J=2 

BOL!REPEAT 

SI!IF(DR.GE.RS(J)) GOTO SAME 
B2!AC=(ACOS(DR/RS(J)))*%(PI7180) 
33!SX (J) =SA(J)=CRS (J) &® 2% AC-DR* (ORS (J) **®2-DRe*2) we, 5) 
B84!GOTO SKIP 

SS!#SAME CONTINUE 
B6!SX(J)=SACJ) 

ST!EH#¥SKIP CONTINUE 
B88!DUCJ)=SX(J) 

89 !jJoJ+l 

9O!UNTIL(J.EQ.2)) 

91!GOSUB INTG 

92! DM=ANS/35 
O3'TFCDM.GE.LQO.AND.OM.LE.HI) GOTO MOM 
O4!IF(DM.GT.HI) GOT MINUS 
95!DR=DR+.1*DR 

96!GOTO LCBLP 

OT!#MINUS CONTINUE 
Q93!DR=NDR-.1*DR 

99!GOTO LCBLP 

}OO!#MOM CONTINUE 

lO1!J=2Q 

102!REPEAT 
TOSFpUCI)=E5CI7*SXxX5) 

104! J=J+!] 

IOS'tUNTIL(J.£E0O.21) 

106'GOSUB INTG 
107!LCBNSC=ANS/(DM*35) 
1O8B!DX=CINTV-2)75 

1OOVEX(1)=.5 

1IO!LxX(2)=! 

Lies) =1.5 

112!LX(4) =2 

PIZB!REPEAT LOOP2:J=5.1.¢5.GT.3) 
P14!EX(J) e060 J—-4)*NX+2 

11S!#LooP2] CONTINUE 

DVS USEPEAT EXTRA? l=1.1.0J3.5T.8) 





ey SX LELXC J) 

118!GOSUB RAD 

119!RX(J)=RANS 

120! EXTRA 

121!GOTO GOUT 

122!#RAD CONTINUE 

Peg (XE. LE.LF) GOTO FRONT 

124! 1 F(XL.~GT.LF.AND.XL.LE.LL) RANS=R 
Peore(XE.G)T.LL) GOTO AFT 

126!GOTO LOOP4 

J27!#FRONT CONTINUE 

128! XX=LF-XL 

LOG! RANS=Re ( (1-0 XX/LF) ®®NF) e® (1 /NE)) 
130!GOTO LOOP4 

13)!#AFT CONTINUE 

132! XX=XL-LL 

I33!RANS=R*( 1-(XX/LA) **NA) 

134! #4LOOP4 RINSUB 

I35!#INTG <#SIMPSON“S RULE INTEGRATION 
136!0DD=0 

I37!REPEAT ODDL?:J=2,2.(J.GT.20) 

1338! QODD=ONN+4*DUl J) 

I 39!#0ODDL CONTINUE 

140! EVEN=0 

141!REPEAT EVENL?:J=3.2,(J.GT.19) 

142! EVEN=EVEN+2*9U (J) 

143! #EVENL CONTINUE 
144!ANS=(INTV/3) # (DUC 1) 4ODN4EVEN*+DU( 21 )) 
145!RTNSUB 

146! 4#GOUT 

147!LD=LOA/D 

148!PRNT 

149 !PRNT 

ISO!PRNT ENVELOPE DISPLACEMENT: [(DENV] 
IS|'!PRNT LENGTH OVERALL: (LOA) 

152!PRNT OVERALL PRISMATIC: [(CP] 
153!PRNT DIAMETER: [(D} 

}54!PRNT LENGTH/NIAMETER: (LD) 

I55!PRNT MIDBODY LENGTH: [(LMB) 

IS6!PRNT FAWN PRISMATIC: [(CPF} 

157!PRNT AFT PRISMATIC: [(CPA]} 

158!PRNT ENTRANCE? (LF] 

ISS!PRNT RUN: [LA] 

160! PRNT LCB SUBMERGED: [(LCBSUB} 
161!PRNT LCB @ NSC : (LCBNSC] 
162!PRNT DRAFT @ NSC : [(R#DR} 

163!PRNT SATISFACTORY (S) OR RECALCULATE (R)? 
164!REAND &DECID 

16O5!IF (&DECID.EQ@."R") GOTO BRANCH 
166!PRNT PLEASE WAIT WHILE OUTPUT AND PASS FILES ARE WRITTEN. THE OUTPUT 
1O7!PRNT FILE FOR GEOMETRY SUMMARY DATA WILL BE “MIT. GEQMOUT". A FILE OF 
168!PRNT OFFSETS WILL BE iN "MIT.OFFOUT". 
L69!QOPENW 3,"MIT.OFFOUT" 

170!3Q="LENGTH WRT FP OFFSET" 
171!GOSUB W 

172!LS(1)=0 

[/3'RS C120 

174!&8X=L5¢1) 

1/5'!8C=RS C1) 

1749R5=" " 





ere 


177!&Q=8X4+8S48C 

178!WRITEF 3,80 

17O!REPEAT LOOPS: J=1.1.(J.GT.8) 
1380 !'!&X=SLX( J) 

181!&C=RX(J) 

182 !8Q=8X+8S+8C 

183!GOSUB W 

184! Jase 

1385!¥€#LOOP5 CONTINUE 

IS6!REPEAT LOOP6:J=2.1.(J.GT.21) 
I3d7!&X=LS (J) 

1388 !&C=RS( J) 

189 !8Q=8&X+8S+8&C 

190!GQOSUB W 

ISU! #LO0OP6 CONTINUE 

192!GOTO FILE 

1I93!4¢W CONTINUE 

1I94!IWRITEF 3.80 

I95!RTNSUB 

1I96!#FILE CONTINUE 

IS 7!OPENW 4,"MIT.PASSHP" 
198!&X=CP 

IOO'WRITEF 4,8X 

200! &X2LO0A 

2OIIWRITEF 4,8&X 

202!8X=D 

Q2O3B!ANRITEF 4,8X 

204 !&X=SF 

2OS!IWRITEF 4,8X 

2U6IWRITEF 4,8NAME 

207!0PENW 5.“°MIT.GEQMQUT" 
2O8!RX="GEQMETRY OUTPUT FOR "“+&8NAME+" 34 
2OG!IWRITEF 5,AX 

CUO ARITer 5,85 

2lI!WRITEF 5,85 

212!'&8P=LOA 

213!8X="LENGTH OVERALL: "+8&P 
2I14!WRITEF 5,&X 

215!&8P=LMB 

216!8X="MID-BODY LENGTH: “"+&P 
Jl Weiler 5.,8X 

218!&P=CP 

2ISI&X=40QVERALL PRISMATICs #+8P 
220!WRITEF 5.&X 

221) !&8P=LN 
222!8X="LENGTH/DIAMETER: "+4+8P 
223!WRITEF 5,8&X 

224!8&P=LCBSUB 

225!1S8X="SUBMERGEN LCB (AFT OF FP): “48P 
P2O!WRITEF 5,8&X 

2271 TT=R+DR 

223 !8P=TT 

229!8X="DRAFT @ NORMAL SUQF. COND.? “"+&P 
230!WRITEF 5,&X 

231!8P=LCBNSC 

232!'!&X="_LCB @ NORMAL SURF. CONN. : "+&P 
233 WRITEF 5.&X 

234!WRITEF 5.AS 

235!RX="RECAP OF INPUT DATA3" 
23O!WRITEF 5 .AX 





23 7!WRITEF 5.85 
238!8P=DENV 
239!8X="ENVELOPE DISPLACEMENT :24+&P 
240!WRITEF 5,8&X 
24)!&8P=LF 

242! &X="ENTRANCE® "+4&P 
243!WRITEF 5,8&X 
244!8P=LA 

245!&8X="RUN 2 44+8p 
246!WRITEF 5,8&X 
247!8P=N 
243!RX="HTAMETER: “+8&P 
249!WRITEF 5.8&X 
250!&P=CPF 
251!&X="FWD PRISMATIC: “+&P 
2S2°WRITEF 5.&X 
253!&8P=CPA 

254!8X="AFT PRISMATIC: "+&P 
255!WRITEF 5.8&X 

256! 0PENW 6,*MIT.LCB” 
257!&X=LCBNSC 
OOS!WwRITEF 6.&X 
259!&X=LCBSUB 
260!NRITEF 6,8X 

261 !8&X=R 

262!WRITEF 6.8X 
263!8X=L0A 

2O4!IWRITEF 6,8X 
265!0PENW 7,"MIT.PASSPHI" 
266! 8&X=) 

267!WRITEF 7,8&X 

268! &X=SL0A 

269!WRITEF 7.8&X 
270!&X=LF 

2T\!WRITEF /7.&X 
272!&X=LA 

273!WRITEF 7.8X 
274!&X=LCBSUB 
CISeNRITEP 7.8&X 
2/6'!DSUB=DENV-DFF 
277!&X=DSUB 

27B8!WRITEF 7,.8&X 

279 '!&X=DENV 

230!WRITEF /,&X 
2381!&X=NF 

262!WRITEF 7,8&X 
2B3!A&X=NA 

O2S4!WRITEF 7,&X 

235! #DRAW 

236!X(1)=0 

287!Y(1)=0 

288!]=2 

289 !REPEAT 
Q290!XCT)=LXC 1-1) 
291!YC1) =RX( 1-1) 

292! [=I +1 
OO3!UNTILCI.GT.9) 

294! LFINTH=(LFE-CINTV+4+1 01723 
295! XL=INTV+! 

296! T=10 





2OT!REPEAT 

293! X(1)=XL 

299 !GOSUB RAD 
300!Y(I)=RANS 
301!XL=XL+LFINT 

302! =I +1 

Bose UNIILEC!.GT.31) 
304!LINT=(LOA|LF)7101 
305!XL=LF*¢LINT 

306! =32 

307! REPEAT 
3O08!X(1)=XL 
309!GOSUB RAD 
310!YC(1)=RANS 
se xe =Xl+LI NT 
312!'l=Il+1 
313!UNTILC(CI.GT.131) 
314!X(132)=L0A 
315!Y(132)=0 
316!0PENW 8."MIT.POINITS" 
317!J=1 

SIS!REPEAT 

319 !8X=XK(J) 
320!8&Y=Y( J) 
321!WRITEF 8,8&X 
B322!WRITEF 8,&yY 

323! J=J+1 
324!UNTIL(J.GT.132) 
325!RUN NEW MIT.SPEED 
326!END 


MIT.&BCD.SPEED 
4-11-84 13:20:34 FUTIL 6.18 


1!IDIM HP(40),KW(40) 

2'READCENTER HOTEL AND COMBAT SYSTEM LOADS IN KW: )HL 
3!READCENTER PROPULSIVE COEFFICIENT: )PC 
4!QOPENR 1,°MIT.PASSHP# 

SERGADF | ,&éCP 

6!READF 1,8&LO0OA 

T!READF 1.80 

S!READF 1.4&SF 

9!READF 1,&NAME 

IO!CP=&CP 

II tLOAS&LOA 

12!D=8&D 

13!SF=8SF 

14! DA=|.09065E=3 *LOAsD4+11.25 

I5!'J=l 

I6!REPEAT 

I7!RE=( Je}. 689"L0A)/1.27908E-5 
I8!CF=.075/(LG(RE)—2) «#2 
IS'!CR=CFe(1.5e(D/ LOA) *#1.547%(0/L0A) **#3)+.002%(CP-.6) 
20! CT=CF+CR+. 00025 
2I'HP(J)=C1/PC)#. OO872%J x23" (CT#SF+DA) 
22! KW(J)=. 745 7#HP (J)4HL 
Zoe w= | 
24!UNTIL(J.GT. 40) 
25!PRNT 
26!PRNT 
27! PRNT 
28!PRNT SPEED AND POWER QUTPUT FOR [&NAME]) 
29! PRNT 
30!PRNT NON PROPULSION LOADS (KW): (HL] 
31! PRNT PROPULSIVE COEFFICIENT EPC} 
32!PRNT 

33!PRNT HIT <RETURN> TQ SEE OUTPUT VALUES... 
34!READ &CONT 

35!PRNT 
36!PRNT SPEED (KT) SHP TOTAL KW 
37!PRNT 
38!J=I 

39 !REPEAT 

40!PRNT [J] (HP (J) ) (KW(J)) 
41! JjJoJ+) 

42!UNTIL(J.GT.20) 

43!PRNT HIT <RETURN> TO CONTINUE... 
44!READ &CQONT 


45!PRNT SPEED (KT) SHP TOTAL KW 
46!PRNT 

4/!REPEAT 

438!PRNT [J] (HP(J)) (KWC(J)) 
49! J=J+| 


SOUCUNTIL(C J) .GT.40) 

51!PRNT PLEASE WAIT WHILE THE OUTPUT FILE IS CREATED. OUTPUT FOR THIS 
S52!PRNT MODULE WILL BE IN FILE "MIT.POWER'". 

SSIQPENW 2, "MIT. POWER* 

54!'&8X="SPEED AND POWER RESULTS FOR "“+&NAMES3 " 

SS5!WRITEF 2,8&X 

56!&PADI =" “ 


5 7!&PAD2=* = 
SS8!WwRITEF 2.&PAD! 
59! 8HL=HL 
60!&PC=PC 


611 &Xs"*NON=]PROPULSION LOADS (KW) 8 #48&HL 
62!WRITEF 2,8X 

63!&X2"PROPULSIVE COEFFICIENT s #+&PC 
64!WRITEF 2,4X 

OS!WRITEF 2,&PADIi 

66!&XS"SPEED (KT) SHP TOTAL KW" 
6O7!IWRITEF 2,&PADI 

68!WRITEF 2,8&X 

69!WRITEF 2.&PAD1 

70! Jel 

T1!REPEAT 

T2!&JaJ 

73!&X=2HP (J) 

T4!QHP=8X(1,8) 

T5!&Y=Skw (J) 

T6!&KW=8Y( 1,8) 

TT! &OUT=" *~+&J+&8PADI +" *+R8HP+&PAD2+&KW 
T8!WRITEF 2,80UT 

719! JaJed 

SO!UNTIL(J.GT.40) 

81!RUN NEW MIT.ENV 

B2!END 


=O) = 


MIT 268CD. BALANCE 
4-27-84 10:47:27 FUTIL 6.21 


ibiwen (oO) seec( 16) .VC(16) .LMC16) .VMC16) .VYC16) 
2!QOPENR 1,2°MIT.BALPASS* 

3!J=| 

4!REPEAT 

S!READF 1,&X 

6!w( J) =&X 

7! jeje 

B8!UNTIL(J.GT.16) 

O!OPENR 2,."MIT.LCB" 

IO!READF 2,&X 


11 'LCBNSC=&X 
I2!READF 2,&X 
13!LCBSUB=&X 
I14!READF 2,8X 
I5!'R=&X 
16!READF 2,8X 
1 7!LOA=&X 


I8!PRNT THE INPUT DATA FOR THIS MODULE AND THE POLYGON MODULE WILL BE 
19! PRNT PROMPTED IN A NEW FORMAT. WHEN THE PROGRAM PROMPTS FOR AN INPUT 
20!PRNT VALUE, THE CURRENT VALUE WILL BE DISPLAYED TQ THE RIGHT OF THE 
21!PRNT INPUT PROMPT MESSAGE. THIS WILL APPEAR AS : 

22!PRNT 

23!PRNT “INPUT PROMPT: = CURRENT VALUE’% 

24! PRNT 

25!PRNT THIS WILL FACILITATE CHANGES TQ INPUT DATA IF MORE THAN ONE ITER- 
26!PRNT ATION THROUGH THE MODULE(S) IS REQUIRED. 

27! 45TART CONTINUE 

23!LT=0 

29! VT=0 

30!J=l 

31!REPEAT 

32!PRNT 

33!PRNT GROUP (J] LCG: 

34!READ(VALUE?)LC(J) 

35!PRNT GROUP (JJ VCG: 

36!READ( VALUE?) VC(J) 

37!LM(JISW( J) «Lo(J) 

3B8!VM(J)=SW(J)*VC(J) 

39! LT=LT+LM (J) 

40!VI=VI+VM(J) 

41 !J=J+l 

42!UNTIL(J.GT. 7) 

43!LC( 3) =LT/W(8) 

44!VC(3d)=VT/W (8) 

45!LM(8)=LT 

46!VM(3)=VT 

47!PRNT 

48!READ(LCG OF VARIABLE LOAD?)LC(11) 

49!READ(VCG OF VARIABLE LOAND?)VC(11) 

50!PRNT 

SIUILMC LIV SWNC 11) *LCOCI11) 

S2!IVM(C TI SWC HET) * VCC 11) 

53!PRNT TOTAL LEAD IS (N(9)) TONS. 

54° PRN 

SS!REAND(CAMOUNT OF LEAD FOR MARGIN IN TONS?) ML 

56!SL=W(9 )=-ML 


2) On 


57TI!LSL=(W0 12) *LCBNSC-LM(8)=-LM(0.11 )-ML*®(LOAZ2 ))/SL 
58!LC(9) =(SL*®LSL#ML*®(LOAZ2))/N(9) 
S9O!PRNT REQUIRED STABILITY LEAD LCG IS [(LSL]. 
6O!REANCIS THIS FEASIBLE <Y/N>? )&R 
6I1!IF(&8R.EQ."N") GOTO START 
62!LM(9)=WN(9)*LC(9) 
63!VM(9)=SL*6+ML*R 
64!VC(9)2VM(9)/WN(9) 

65!LM(10) =LM(8)+LM(9) 
66!VM(10)=VM(8)+4+VM(9) 
67!LC( 10) =LM(10) WC 10) 

68! VC(10)=VMC10)7N0 10) 
69!LM(12)=LM(10)+LM( 11) 
TO!VM( 12) =VM( 10) 4VM(11) 

71!VC(13)=R 

7T2!1VM (013) =R*¥W013) 

73!LCC13 =(WO014) *LCBSUB=-LM(12))/N(013) 
74! PRNT 

75!PRNT REQUIRED MBT LCG IS [LC(1I3)]. 
T6!READCIS THIS FEASIBLE <Y/N>? )&R 
TT! TFCA&R.LEQ."N") GOTO START 
73!LM(13)=W013)*LC(13) 
79!'LM(14)=LM(12)4LM(13) 
BO!LC(14)=LM(14)7W014) 
81!VM(14)=VM(12)4+VM(13) 
B2!'VC( 14) 2VM014)/N014) 

83! BG=R-VC(14) 

B84!PRNT 

85!PRNT SUBMERGED STABILITY (BG) IS (BG) FT. 
B6!READCIS THIS ACCEPTA3LE <Y/N>? )&R 
S7TIIFCAR. EQ. "N4) GOTO START 

88!PRNT THE OUTPUT FOR THIS MODULE WILL BE IN FILE “MIT. BALOUT". 
B89! J=] 

SO! REPEAT 

91'VYCJ)SVCCJ)—R 

92! JoJ+l 

O3!UNTIL(J.GT.14) 

O94!'Q0PENW 3.°MIT.BALOUT" 
95!&8X="RESULTS FROM BALANCE MODULE" 
O6!WRITEF 3,8X 

97!'&P),=" at 

98!&Pas" % 

99!'!8P3=s" «I 

IOO'IWRITEF 3.4P! 
101!&X="*GROUP WEI GHT LCG Veo" 
IO2ZIWRITEF 3.&X 
IO3!WRITEF 3,8P1 

104! J=] 

105! REPEAT 

106!&J=sJ 

1O7!t&X=en( J) 

108 !8w=&X (1,8) 

1O9'&X=LC(J) 

110!&LC=&X(1 .6) 

111!a@x=sve(J) 

112!'&VC=8X(1 66) 

11 3!8OUT=" 4448 J++ N+BN+8P2+ALC+A&P3+8VC 
LI4¢WRITEF 3 ,8&0UT 

115! jJaJel 

LrotUNta oC S.CT.) 





[ivi Jaa 
118!&G=" Al i 
119!GOSUB [0 

1 20!J=9 

VIS G=" LEAt * 
}22!GOSUB [0 
l2g5J=10 
124!&Gs" A a 
125!GOSUB IQ 


126!J=11 
127%&G=" V LD " 
128!GOSUB I0 
129!J=13 
130!&G=" MBT * 
131!GOSUB IQ 


1} 32!&NSC=W( 12) 

133!&SUB=N (14) 

134!&LCBN=LCBNSC 

135!&8LCBS=LC BSUB 

136!&BG=BG 

lS TIWRITEFP 3,&P1 

138!&SL=SL 

139 !&8ML=ML 

140! &8QUT="MARGIN LEAD (TONS): “+#&ML4¢" @VCG=D/2, LCG=LOA/2" 
IS41!WRITEF 3,8QUT 

142!8OUTS"STABILITY LEAD (TONS): *%4&5L44* @VCG= 6 FI" 
143!WRITEF 3,80UT 

144'&8LSL=LSL 

145!&OUT="STABILITY LEAD REQUIRED LCG: "“+&LSL 
E4o ARI T EF 32 ,£0UT 

IS4/!WRITEF 3,8P1 

148!R8OUTS“NSC (CTONS) 2 "“#&NSC+44 @LCB: “+&LCBN 
}49!WRITEF 3.8QUT 

ISO!WRITEF 3,.&P1 

151 '!&8OUTH="SUBM (TONS): “4+&8SUB4" @LCBs: “+&LCBS 
ISQ2IWRITEF 3 &0OUT 

IS3'WRITEF 3,4&Pl 

154'&OUT="SUBMERGED STABILITY (BG): "+#&BGe" FT4 
loosen Gree 3,5 0UT 

I56!0PENW 4,."MIT.PASSPOLY" 

157!&8X=LCBSUB 

ISS!wRITEF 4,&X 

159 !&8X=LC( 11) 

16O0!WRITEF 4,8X 

161!&XsLC(13) 

I62'WRITEF 4,8&X 

163!GOTO PLOT 

164!#1Q CONTINUE 

165!&X=w( J) 

166!8W=8X(1,8) 

lO 7S X=LC (3) 

16O3d!&LC=AX( 1.6) 

16GVEx—VCC J) 

17OTEVEH=4AXC1 26) 

171 $8OUT=&8G+A8W+8P2+&8LC+&P34+8VC 

1} /2!WRITEF 3,a&0UT 

173!RTNSUB 

174!#PLOT 

1¢S5'SELECT MODE MODEL 

1 /6!SELECT LAYER 3 





P77! J=I 

| 78!REPEAT 

179%8&J=J 

ISO'INSERT POINT? XELC(J)IY(VY(J)1<CR> 

toue *P=LeEC J) +1 

[oeeNoehn Text (81) HGT 1.52? X(XPIJY(VY (J) 1]<CR> 
I83°J=J+1 

IS84!UNTIL(J.GT. 7) 

IB8S!tINSERT POINT? X(LC(9) J Y{VY(9) ]<CR> 
IB86!XP=LC(9) +1 

I87!INSERT TEXT “LEANDY HGT 1.53X{XPIJY{VY(9)]<CR> 
I8s!INSERT POINT®X(LC( 13) JYCVY (13) 1<CR> 
I89!XP=LC(13)+1 

[ISeMINSERT TEXDeeMb)’ HOT 1.58X(XPJY(VY(13)1<CR> 
I9S1!RUN NEW MIT.POLY 

IS2!END 


MITZSBCDSPOLY 
5- 77-84 19214:02 FUTIL 6.21 


ine meron ClO), YV (10) .XP(10) .YP(10) .XTC10) .YT(10) .M(3d) ,Y(8) 


2'QOPENR 
3!J=1 
4!REPEAT 
S!IREADF 1,&X 
6!W(J)=8X 
(ljej+i 

Ss UNtrCJ.GT. 16) 
OIVLN=EWC 11) 
}O!4BT=W (13) 

Peo SsuB=wCi4) 


ee he BALPASS* 


[200PENn ee. Mi 1.PASSPOLY" 


I3!READF 2,&X 
14!LSUB=8X 
IS5!READF 2.&X 
1O!LVLD=&X 
| /!READF 2,&X 
IS !LMBT=&X 


19! #START CONTINUE 


20!BAL=VLD 


Ai eon DOmAL VARIABLE LOAD 15 (VLD) TONS. 


22! PHNT 
23'READCWEIGHT OF 
24! READ (LCG?) LCE 
25!BAL=BAL=-CE 
26'O0SUB STAT 
27!REAN(CNEIGHT OF 
25!READ(LCG?) LMC 
29!BAL=BAL~MC 
30!GOSUB STAT 
31!READCWEIGHT OF 
32!REAI( LCG?) LST 
33 !BAL=BAL-ST 
34!GOSUB STAT 
35!READ(CWEIGHT OF 
36!READCLCG?) LLO 
3 /!RBAL=BAL-LO 
38!GOSUBM@STAT 

39 !READ(WEIGHT OF 
40! READ (LCG?) LCFO 
41!BAL=BAL-CF() 
42!GOSUB STAT 
43!REAN(WEIGHT OF 
44!READ( LCG?) LPS 
45!RAL=BAL=—PS 
46!GOSUB STAT 
47!READCWEIGHT OF 
S48! REA) (LCG?) L02 
49 !BAL=BAL-02 

DUE GOSUS SPAT 
SItRFANCWEIGHT OF 
52 HEADIN LCG?)LTF 
SS3YBALHEBAL=fF 
S4GOsusS. STAT 
SOS!'Ye AN (WEIGHT OF 
54! aEANCI.CG?) LTA 


SHEWeAND ErFECTS?)CE 


OR COMPENSATING WATER?) MC 


SLBM%5 


SANITARY TANKS AND WATER?)ST 


LUBE OIL IN SUMPS?)LO 


SHIELD TANK>?)CFO 


<IE, 


RIXER CLEAN F. 


Q. 


PROVISIONS AND STORES?)Ps5 


REVITALIZATION OXYGEN?) 02 


TORPEDOS IN FORWARD ROOM?) TF 


TORPEDOS IN AFT ROOM?)TA 





07 = 


57!BAL=BAL-TA 

SSo!'GOSUB STAT 

SO!READCWEIGHT OF TACTICAL MISSILES FORWARD?) MF 

OU!LREAN( LCG?) LMF 

61!BAL=BAL-MF 

62!GOSUB STAT 

63'REAND(WEIGHT OF TACTICAL MISSILES AFT?)MA 

64'READ(LCG?)LMA 

o5!BAL=8AL—=MA 

66!GOSUB STAT 

O7!TMSTF+eTA+MF+a 

6G6!LTM=(TF*eLTFeT A®*LTA+MFexLMF+MAeLMA)/STM 

OD! READCWEIGHT OF WATER FOR WRT TANKS?) WRT 

/O!READ (LCG?) LWRT 

71 !BAL=BAL=WRT 

f2'GOSUB STAT 

7T3'READCWEIGHT OF RESERVE ELECTROLYTE?) EL 

74!READ (LCG?) LEL 

75!BAL=BAL-EL 

fo'!GOSUB STAT 

TT! READCWEIGHT OF FRESH WATER <POTTABLE.FEED.BATTERY RESERVE>?)FW 
73 !READCLCG?) LEW 

/9'BAL=BAL-FW 

sO!GQOSUB STAT 

SI'READCWEIGHT OF RESERVE LUBE & HYDRAULIC OIL?) RLO 

32!READ( LCG?) LRLO 

33!BAL=BAL-RLO 

384!GOSUB STAT 

85!PRNT THE LAST VARIABLE LOAN ITEM TO BE CONSINERED IS FUEL OIL. TO ENSURE 
S6!PRNT A MARGIN OF SAFETY AT THE BASE OF THE POLYGON, A VARIABLE FUEL OIL 
S7TIPRNT (VFQ) TANK IS AUTOMATICALLY SPECIFIED. THE REQUIRED VFO CAPACITY IS 
83! PRNT COMPUTED USING YOUR INPUT OF THE FRACTION OF TOTAL FUEL TO 
SO!PRNT PLACE INTO THE VFO TANK. THE MINIMUM SUGGESTED FRACTION IS 
QUIPRNT .23 (THE RATIO DIFFERENTIAL OF FUEL AND WATER DENSITIES). 
QI'PRNT IT IS HIGHLY RECOMMENDED THAT AN AMITIONAL MARGIN OF AT 
Y2'IPRNT LEAST .05 BE SPECIFIEl) TO ALLOW A SAFETY MARGIN AT THE BASE 
Q3!PRNT OF THE POLYGON. 

Q4!PRNT 

Q5!REFADC FRACTION OF FUEL TO Be PLACED IN VFO?) VER 

YO!TFO=BAL 

QO 7! VEO2VER* TFO 

YQ3!FBT=TFO-VFQ 

99! FL=ECE+MC+$T/34L0+CFO 
100! LFL=(CE*LCE+#4C*LMC+ST*LST #L0* LLO*+CFO*LCFO) /ZFL 
1OlT!PRNT FBT CAPACITY (TONS): [FBT] 

IO2!PRNT VFO CAPACITY (TONS): (VEFQ} 
103! PRNT 
LO4!*READCAMOUNT OF FBT <TONS> IN FORWARD TANKS?) FBTF 
}O5!FBTA=FRT-FBIF 
}O6!PRNE 

1O7!READCFWD FRT LCG?)LFSATF 

IOS! REAICAFT FBT LCG?)LEBTA 

LOO! LEFRT=(FRTF*eLFBTF+FBTA*LEBT 4) /FBT 

ITO!REANCVEQ LCG?) LVFO 

PITT EPRNE 

11 2'PRNT 

11 3!PRNT 

11 4!PRNT 

}lis!PRNT VARIABLE LOAI) RECAP: 

1}16!PRNT Ife WE 1GHT EGG 


=\oje\a 


11 7! PRNT 

[Tease RNieeCREW & EFFECTS [CE] Gree 
119! PRNT SLBM7S (MC } (LMC } 
}2Q0!PRNT SANITARY TANKS (st) PES 
121)!PRNT SUMP LUBE OIL (LO) CLO 
P2Se ee GEEAN FE. . (FIXED) (CEO) GLCEO! 
123!PRNT ***TQTAL FIXED*** bebe) GREE 
124!PRNT PROV. & STORES Ons) LEVS) 
I25!PRNT REVITALIZATION 02 (02) (LO2Q] 
126!PRNT TORPEDOS FWD ROOM aes) Cethy 
}27!PRNT TORPEDOS AFT ROOM (TA] Ghee 
P2SUPRNG TCT. MISSILES FwO (MF } (LMF } 
}29!PRNT TACT. MISSILES AFT (MA] ( LMA} 
130!PRNT WRT TANKS (WRT (LWRT] 
PS) SPRNT RESERME ELECTROLYTE (EL!) Piers 
Pa2-rniny TORAL, FRESH WATER (FW) (LFW) 
WSS "PANT RESERVE LUBE OIL (RLO} (LRLO] 
134!PRNT FUEL BALLAST FWh CEBEP) CEEB EP 
135°PRNT FUEL BALLAST AFT (FBTA] (LFSRA] 
I36!PRNT VARIABLE FUEL OIL (VFO) CLVEO) 


I37!READ(CARE THESE VALUES ACCEPTABLE <Y/N>? )8&R 
I3BB!IFCARLEQ."NO") GOTO START 
139!PRNT 

140! AFL=LSUB-LFL 
141! APS=SLSUB=-LPS 
142!AQ2=LSUB-LO2 
143!ATM=LSUB-LTM 
144! AWRT=LSU B-LWRT 
145 ¥AEL=LSUB<-LEL 
146! AFW=LSUB=-LFEW 
147! ARLO=LSUB<-LRLO 
140! AF BT=LSUB-LFBT 
149TAVFQO=LSUB<-LVFO 
1SO!ATFSLSUB-LTF 

151 !ATA=LSUB-LTA 
152! AMF=LSUB=LMF 
153! AMA=LSUB=-LMA 
154!AFBTF=LSUB=-LFBTF 
155! AF BTA=LSUB=LFBTA 
156!MFL=FL*AFL 
157!MPS=PS*APS 
153!M0 2=02eA02 

159! MTM=TM®ATM 

150! MWRT=WRT *AWRT 
161 !MELS=EL#AEL 
162!MFW=FWeAFW 
163!MRLO=RLO#®ARLO 
164!MFRT=FBT*AFBT 
165!MVEFO=VFOeAVEQ 

| 60! MATF=TFeATFE 
1\O7!MTA=TA*®ATA 

log! MMP=MF RAMP 
|\69!MMA=MARAMA 

1 7U!MPBRTF=FBTFeAFBTF 
| /1'M4FEBTASFBTA*®AFBTA 
172!CVL=.99375 

1 (Sev H=120093 7/5 
174!GOTO NSTM 

1 /5!e@STAT CONTINUE 

| /6!PRNT 


1 77'!PRNT BALANCE OF VARIABLE LOAD ACCOUNT: [(BAL) 

173 !PRNT 

1/79!RTNSUB 

180! #NSTM CONTINUE 

Isl! LOSUB=H=N (14) =-W (10) 

1B2!#HEAVI CONTINUE 

183!CW=43/35 
IS4!HI=FL+TM+.5*PS+EL4+wRT +. S* RLO+.5*FW+CWe FBT 

185! MHI =MFL4+MTM+. S*MPS*MEL+MWRT +. 5 ®MRLO+, S*® MEW +CNeMEBT 
136! VBH1=LDSUB-—MBT-H1 

IS7!AMBT=LSUB=-LM8T 

| 88!MMBT=MBT*AMBT 

139 !MVBHI =—MMBT-MH1 

190! XV (1) =MVBHI*CVL 

191! YVC1)=VBH1 *CVL’ 

OP! #HEAV? CONTINUE 

193'H2=FL+TM+PS+EL 4¢WwRT+RLO+FW+4CW * FBT 
194!MH2=MFL4¢MTN+MPS4MEL*+4NRT*MRL OMEN +CW eM PBT 

195! VBH2=LDSUB-—MBT-H2 

196!MVBH2=-{MMBT -—MH2 

197! XV (2) =MVBH2*CVL 

IS38!YV(2) =VBH2*CVL 

199! @LITE} 

200! LISFL+. 75*PS4WRT +. 75*RLO+. (5 *FW+FBT+VEQ) 

201 !ML1=MFL+,. 75*MPS4#MWRT +. /5*MRLO+. /5*MPN+MFBT MVE) 

202! VBLI=LDSUB-MBT-LI 

203!MVBL1=—{MMBT-ML1 

204! XV(3)=MVBLI *CVH 

2U5!YV( 3) =VBLI*CVH 

206! #LITE2 

207!L2=FL4¢WNRT+RLO#VEO+Cwe FBT 
2U38!ML2=MFEL+MWRT+MRLOSMV FO +CWeMEBT 

209! VBL2=LDSUB=MBT -L2 

210!MVBLO=-]MMBT=-ML2 

211!XV(4)=MVBL2*CVH 

212!YV(4) =VBL2*CVH 

213! #HF 1 
214!SHFISFL+MF+TF+. 75*PS+EL+WRT +. 75 *RLO+. 75% FW4+VE()*+F RTA+CW®FBTF 
215!MHF1=MFL4MTF+MMF +. /5*MPS+#MEL+MWRT +. (5 *@ MRL. (5S *MEW eM VEFO4MFBIA+CW eM FBTE 
216! VBHF1=LOSUB-MBT-HF 1 

217!MVBHF 1 =-MMBT=MHF 1 

218!XV(5)=MVBHF 1 

219!'YV (5) =VBHFI 

220! #HF2 
P2\SHFQ2=FL+MF4TF+. 5*P> +EL+WRT+. 25*RLO+. Se FW4+VEO+FRTAt+CWeFE BTF 
222!MHF2=MFL+MTF+MMF+, 15 *MPS+MEL +MWRT +. 25 *MRL()+. Se MEW4+MVE04+MEBIA+CW eM FBTF 
223! VBHF2=LDSURB-MBT-HF2 

224 !MVBHF2=~-MMBI=-MHF2 

P25!XV(6) =MVBHFO 

226'YV (6) =VBHF2 

227! #HAa 

P23 !HA =FL+MA+TA+. O*PS4ANRT+. 75 e#RLO+. Se FW 4VFO+FBTF+CWeF BTA 
229! MHA=MFL4+MM AMT A+. 5 ®MPS+MWRT +, 75 e@MRLI+,. 5S ®MEW EMV IF04+MEBTE+CN®MEBTA 
230! VBHA=LIDSU3=MBT-HA 

231!MVBHA=—MMBT=MHA 

232!XV07)=MVBHA 

233!YV(/)=VBHA 

234!8N 

P35!N=HD 

236!MN=MH2Q 


—!BQ—- 


23/7!'V3N=LDSUB=MBT=N 

238!MVBN=—MMBT—MN 

239! XV(8) =MVBN 

240!YV(3)=VBN 

241!#CONDM CONTINUE 

242!CM=H2-CW*FBT+FBT 

243! MCM=MH2—-CW*eMFBT+MFBT 

244! VBCM=LDSUB=-MBT-CM 

245!MVBCM=={MMBT-MCM 

246!XV(9)=MVBCM 

247!YV(Y9) =VBCM 

24d! <#xxxeexxeazeeee SET UP LOCATION OF TRIM & AUX TANKS «etek keke eka 
249 !PRNT 

250!PRNT YOUR ACTUAL NESIGN MAY HAVE AS MANY TRIM AND AUX TANKS AS YOU DEEM 
251!PRNT APPROPRIATE FOR YOUR PARTICULAR ARRANGEMENT SCHEME. TQ SIMPLIFY THE 
252'PRNT CALCULATION PROCESS FOR THE EQUILIBRIUM POLYGON, YOU ARE ASKED To 
253!PRNT ENTER AGGREGATE CAPACITY ANN LCG FOR THE FORWARN TRIM GROUP. THE 
254!PRNT AFTER TRIM GROUP. AND THE AUXILIARY GROUP TANKS. 
255!PRNT 

25964"CG0OT0 SIZE 

257! €# LOOP 

258!EXECV 

259!READ(LCG OF FORWARID TRIM TANK GROUP?)LFT 

260!REAN( LCG OF AFTER TRIM TANK GROUP?)LAT 

261!READ(LCG OF AUX TANK GROUP?)LAUX 

262! AFT=LSUB=-LFT 

263! AAT=LSUB-LAT 

264! AAUX=LSUB~LAUX 

205!RTNSUB 

266! #SIZE CONTINUE 

26/!J=l 

263!REPEAT 

269!TFCXV(J).GT.5000)GOTO LARGE 
270!1TF(XV(J).LT.-S000)GOTO LARGE 
271!IFCYV(J).GT.500)GOTO LARGE 

2/2! JaJel 

27 atting Pot J.GIT.9) 

274!GOTO SMALL 

2/S!#LARGE 

2TO!INEL PAR MIT.HULLOUT 

Q2T7tEXIT PART FILE MIT.HULLOUT 

2738!XS=1000 

2/9!YS=100 

280!ACT PAR MIT.POLYOUT 

231!ACT DRA POLYGON FORM MIT.LARGEPOLY DRAW POLY 

282!SEL CPL LEFT 

233!DEF VIE GRAPH: XI1Y6,XOYO,X22Y17 <CR> 

284!GOTO GRAPH 

295!45MALL CONTINUE 

2uo DEL PAR MIT.HULLOUT 

237!EXIT PART FILE MIT.HULLOUT 

2380! XS=500 

239! YS=50 

29O!DEL PAR MIT.POLYOUT 

29'SACT PAR MIT.POLYOUT 

292!ACT DRA POLYGON FORM MIT.SMALLPOLY DRAAN POLY? 
293!SEL CPL LEFT 

294!NEF VIE GRAPH: XIIYO.XUYO.X22YI1/ <CR> 

295! #GRAPH CONTINUE 

296! J=1 


| 


2OT7T!REPEAT 

2968! XP (J) =XV(J)SXS 

2990! YP(JVAYVIJIZYS 

300!INSERT POINT: XC XP(J)IJY{YP(J)] <CR> 
301!XT( J) =SXP(J)+.2 

B02 !YT(J) =YP(J)-. 1 

303! JeJ+l 

304!UNTIL(J.GT.9) 

BOSVINS TEX H! HG .2: X{XTCI) YC YT(1))] <CR> 

306!INS TEX H2 HG .2: X(XT(2))Y{YT(2)) <CR> 

SOT!INS TEX LI HG .2?: X{XT(3))YCYT(3)) <CR> 

SO8!INS TEX L2 HG .2: X(XT(4))Y{YT(4)] <CR> 

309!INS TEX HF! HG .2: X{XT(5) JY{YT(5)] <CR> 

310!INS TEX HF2 HG .2:2 X{XT(6))JY{[YT(6)] <C> 

3Il!INS TEX HA HG 2.2: X{XTC7)IJYCYT(7))] <CR> 

B1l2!INS TEX N HG .2? X(XT(8)-.5)YLYT(3))] <CR> 

Sis ils TEX M HG .22 X(XT(9O)JYCYT(9)] <CR> 

314!GOSUB LOOP 

315!GOTO PLTNK 

316!4CAP CONTINUE 

317!'EXECV 

3I8!READ( CAPACITY <TONS> OF FORWARD TRIM GROUP?) CFT 
319!'READ(CAPACITY <TONS> OF AUXILIARY GROUP?)CAUX 
320!READ( CAPACITY <TONS> OF AFTER TRIM GROUP?) CAT 
321!RTNSUB 

322! #PLTNK CONTINUE 

323!GOSUB CAP 

324!M(1)=0 

325!M( 2) SAFT*CFT 

326!M(3)=M(2)+AAUX *®CAUX 

327!M(4)=M(3)+AAT*CAT 

328!1M(5)=M(4) -AFT* CFT 

329!M(6)=M(5)-AAUX#®CAUX 

330!M(7)=0 

331!Y(1)=0 

332!Y(2)=CFT 

3331Y(3)=CAUX*+CFT 

334!Y(4)=Y(3)+CAT 

335!Y(5)=Y(4)-CFT 

336!Y(6)=Y(5) -<CAUX 

33 7!Y(/)=0 

338!'J=1 

339 !REPEAT 

340!M(J)=M(J)/XS 

B41!Y(JVSYIJISZYS 

342! J=J+] 

S433UNTIL(J.GI. 7) 

344!INS LIN TAG=12 X(MCIIIJYCYCI)).X€MC2) JYLY(2))]. <CR> 
Beat iNae IN TAG=2: XIM(2) J YCY(2)),X0M(3) )YCY(3))]. <CR> 
346!INS LIN TAG=3: X(M(3B)]YCY(3)).X0M (4) )YCY(4)),. <CR> 
347!INS LIN TAG=4: X{M(4) ])Y(LYC4)] .X0M(05) )Y0Y(5))]. <CR> 
343!INS LIN TAG=5: X(M(5)IJYCY(5)) .X0IM(6) J YCY(6)). <CR> 
349!INS LIN TAG=45: X(M(6)) YL Y(6)]).X€M07) JYCYC7)). <CR> 
350!EXECV 

S5I1!READCIS THIS POLYGON SATISFACTORY <Y/N> 7 YAR 
B52'ITF (AR. EQ. "Y") GOTO FINAL 

353!DEL ENT: TAG 1.TAG 2.TAG 3.TAG 4.TAG 5.fAG 6 <C> 
354!EXECV 

355!REALI)(CHANGE CAPACITIES ONLY <tt>, OR LCG“S & CAPACITIES <2> ? )CHC 
356! 1F(CHC.EQ.2) GOTO BOTH 


as 7tiecedG.cO.1) GOTO PLINK 
358!#BOTH CONTINUE 
359!GOSUB LOOP 
300!GOTO PLTNK 
361'#FINAL CONTINUE 
362!0PENW 3.2°MIT.PASS.POLYDAT" 
363!&X=CE 
364!'GOSUB 0) 
305!&X=LCE 
366!GOSUB 0 
367!&X=MC 
368!GOSUB 0 
369 !&X=LMC 
3/0!GOSUB 0 
371 $&x=ST 
3/2!GOSUB 0 
373!&X=LST 
374!GOSUB 0 
375!&X=LO0 
376!GOSUB 0 
377!&X=LLO 
373!GOSUB 0 
3/9!&X=CFO 
3830!GOSUB 0 
381!&X=LCFO 
382!GOSUB 0 
333!&X=FL 
3384!GOSUB 0 
345 !&X=LFL 
3386!GOSUB 0 
3387!&X=PS 
383!:GOSUB 0 
339 !&X=LPS 
390!GOSUB 0 
391 '&X=02 
392!GOSUB 0 
393 !'a&X=L02 
394!'GOSUB 0 
395!&X=TF 
396!GOSUB 0 

39 /'&X=LTF 
398!GOSUB 0 
399!&X=TA 
400'GOSUB 0 
401!4&X=LTA 
402!GQSUB () 
403!&X=MF 
404!GOSUB 0 
405 !&X=LMF 
406!GOSUB (0) 
407!&X=MA 
408!GOSUB 0 
409!'!&X=LMA 
410'GOSUB () 
411°AX=wWRT 
412!GOSUB 0) 
413! 48X=LWRT 
414!GOSUB 0 
415!@€X=EL 
445!SOSUB 0 





AW eX lee 
4138!GOSUB 0 
AIS! &X=Fn 
420!GOSUB 0 
421!&X=LFW 
422!GOSUB 0 
423'&X=RLO 
424!GOSUB 0 
425! &X=LRLO 
426!GOSUB 0 
427!&X=FBIF 
428!GQOSUB 0 


429!&X=LFBIF 


430!GOSUB 0 
431 !&X=FBTA 
432!GOSUB 0 


433!&X=LFBTA 


434!GOQSUB 0 
435! 8X=VFO 

436!GOSUB 0 
437! &X=LVFO 
4328! GOSUB 0 
439!8&X=CFT 

440!GOSUB 0 
441 !&X=LFT 

442!GQOSUB 0 
443!&X=CAUX 
444!GOSUB 0 
445!&X=LAUX 
446!GOSUB 0 
447'&X=CAT 

448!GOSUB 0 
449 !&X=LAT 

450!'!GOSUB 0 


451!PLOT DOT SCALE -95 
452'EXIT PART F 
453!GOTO END 


454° #0 


455!WRITEF 3.4% 


4560 !RTNSUB 
45 7! #END 


458!RUN NEW MIT.POLYFILE 


459! END 


oe 


—|04— 


Met. 2BCD. POLYPIEE 
5—~ 7-84 19322:03 FUTIL 6.21 


I!OPENR 1,."MIT.PASS.POLYDAT" 
2!QPENW 2,."MIT.OUT.POLY" 


3!PRNT OUTPUT FOR THE POLYGON MODULE WILL BE IN “MIT.OUT.POLY?. 


4¢&T="POLYGON SUMMARY OUTPUT FILE® 


SIMRITErP 2.87 

6't&p=a=" W 

GiWRITEF 2.8P 

BIWRITEF 2,&P 

Q!aT=" ITEM 
TOMNRITEF 2.87 

11!WRITEF 2.8P 

IQ2tREADF 1,&X 

I3!READF 1.8&Y 
14!&T="CREW AND EFFECTS 
Daiwhr ler 2.47 

16!GOSUB RD 

1 7!&8T="SLBM7S OR COMP. WATER 
PoUNR ier, 2,41 

19!GOSUB RD 
20!&T="SANITARY TANKS 
Slater 2. oT 

22!GQOSUB kD 
23!8T="LUBE OIL IN SUMPS 
24!WRITEF 2.8T 
25!GOSUB RD 
eo et Fike) CLEAN FUEL OIL 
SIAN RIleEr 2,41 
28!GOSUB RD 


2928TH" eee TOTAL FIXED ITEMS «*«* 


B30!WRITEF 2.&T 

31!GOSUB RD 
32'&8T="PROVISIONS AND STORES 
SS°WRITEr 2.81 

34!GOSUB RD 

BO AT="REVITALIZATION 02 

BOI WRITER 2 AT 

37!GOSUB RD 

33!&8T="TORPEDOS IN FWD ROOM 
SO!WRITEF 2,8&T 

40!GOSUB RI) 

41'&T="TORPEDOS IN AFT ROOM 
Act WRI ieee 6. 

43!GOSUB RD 

44!&T="TACTICAL MISSILES Fawn 
45!wRITEF 2.4T 

40!GOSUR RD 

4/'&T="TACTICAL MISSILES AFT 
AS!IWRITEF 2.,&T 

49!GOSU8 RI) 

SU'STH"WRT TANKS 

D1 WR ieee 2 AT 

52! GOSUB Rf) 

D3'ATSE"QESERVE ELECTROLYTE 
54°WRITEF 2.4T 

55 "GOSu4 Hl) 

S5O0!8&THa" TOTAL FRESH WATE? 


LCG 


4 RY + a! 


WHE RY > 


WHR You 


He RV >a 


He RY ee 


HHeRY Ht 


Ne RY + at 


WH QV +i 


UHR YV +I! 


He QV + iT) 


we RV +e 


WeRYVHet 


WHR Y +4 


NE RYVeu 


WHRY HU 


NelGaT 


"+ BX 


WH aX 


"+X 


4148 X 


W4&X 


"48 X 


"+ &X 


w+ RX 


"+&X 


"4 RX 


+e eX 


Wee X 


"48 X 


Wea X 


+k X 
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S/tWRITEF 2.aT 

58!GOSUB RD 

59'&8T="RESERVE LUBE OIL He RYVou 
6O!WRITEF 2.aT 

61!GOSUB RD 

62!:&T="FUEL BALLAST FWD H+ RY+e 
So WRITEP 2.41 

64!GOSUB RD 

O5!&T="FUEL BALLAST AFT Ne RVs 
So MR ITEF 2.81 

67!GOSUB RD 

Oo8!&T="VARIABLFE FUEL OIL "+ QYou 
OGIARITEF 2,481 

7O!GOSUB RD 

MVIWRITEF 2,8P 

T2!INRITEF 2.8&P 

(3!&T="VARIABLE BALLAST TANK DATAs" 
(4!WRITEF 2.&8T 


75!&T="" TANK GROUP LCG 
POVARITEF 2.aT 
7T7!&T="FORWARI) TRIM H+RYou 


BoA ITE 2.47 
79 !GOSUB RD 


SO'AT="AUXILIARY ray +e 
Saenkirer 2,47. 

82!GOSUB RD 

SStST="APTER TRIM Bee et 


Ba WRITEr 2,47 
85!GOTO END 
S6!€#RD CONTINUE 
S7!READF 1 .4&X 
S8!READF 1,.8&Y 
89!RTNSUB 

90! #END 

SISEND 


HH ax 


"+8 X 


1+ &X 


"HBX 


CAPACITY. (TONS ?" 


N48X 


"+8 X 


442 X 


207066 
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